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FOREWORD 


This volume of the final report is submitted to the MSA Langley 
Research Center in partial fulfillment of Contract MS 1-10973* This 
contract involves the formulation of a mathematical model for predicting 
the dynamic behavior of a class of rotating flexible space -station con- 
figurations, and the preparation of associated computer programs. 

Volume I contains the development of the required theoretical 
techniques, and Volume II is a computer-program user's manual. 

Dr. Robert W. Fralich of the MSA Langley Research Center is the 
Technical Monitor. Dr. Fred Austin of the Grumman Aerospace Corporation 
is the Project Engineer. 

Mr. Eugene J. Lowe prepared all of the computer programs described 
in this manual and wrote the descriptions of the programs and the method 
of setting up the cards. Dr. Fred Austin wrote the engineering descriptions 
contained in this manual. 

The authors gratefully acknowledge the useful comments provided by 
Dr. R. Fralich and Mr. J. Smedfjeld. Also, Miss B. Durling's help in 
converting the computer program is sincerely appreciated. 
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NOTE ON UNITS 


The program may be used with all consistent sets of units. For example, 
if the input data is in the lb -ft -sec -slug system, the output will also be 
in that set of units; if the input data is in the Newton-meter-sec-kg 
system, the program will express the solutions in that system. 

In the Phase I program only, the user may input the weight and weight 
moment of inertia of each mass and the acceleration of gravity, G. The 
program then divides these weight properties by G to obtain the mass propertie 
If the mass properties are known, the user may input these instead of the 
weight properties; however, in this case, a value of G equal to one must be 
input . 
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1.0 INTRODUCTION 


This user's manual provides the information required to use and modify two 
computer programs. The programs were developed on the IBM 370-165 computer 
located at Grumman Data Systems and then converted for use on the CDC 6600 com- 
puter located at the NASA Langley Research Center. 

The first program, known as Phase I, is used to obtain the vibration modes 
of a structure using the modes of the component modules. The structure may 
consist of up to seventeen modules. The second program, known as Phase II, is 
used to solve, by numerical integration, the equations of motion of a flexible 
two-mass Space Station consisting of a Laboratory, a Counterweight, and a Con- 
necting Structure. The modal characteristics and mass distribution of the 
Laboratory and Counterweight are required input data for the Phase II program. 

These properties may be generated by using the Phase I program or may be obtained 
elsewhere. The output format of the Phase I program is compatible with the input 
of the Phase II program, and the two programs can be run together in one sutmittal. 

The Phase I program requires 70,000 octal core-storage locations, and the 
Phase II program requires 300,000 octal core-storage locations on the CDC 6600 
computer. 
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2.0 PHASE I — MODAL COUPLING OF SUBSTRUCTURES 


2.1 INTRODUCTION 

The Phase I computer program is used to obtain the modes of vibration for a 
structure using the modes and mass properties of each of its component sub- 
structures. Figure 2.1 illustrates, schematically, the most general configuration 
■which may be treated. There is a maximum of seventeen substructures. Also, less 
complex configurations may be considered since any of the substructures may be 
absent from the idealization provided, of course, that the absent substructures 
do not disconnect the configuration. The modules are shown in-plane and 
perpendicular to each other for schematic purposes only. Actually, they may be 
skewed at any orientation in three dimensions. In fact, substructure modes can 
be supplied in a coordinate system which is skewed relative to the total-structure 
coordinate system. The central modules (numbers 1, 11, 14, 16 and 17) are called 
core modules, and the outer modules are called appendages. 

While the total -structure modes to be used as input data for the Phase II 
Program are free-free, the Phase I Program may also be used to obtain cantilever 
modes for the total structure with Module 1 fixed at some point. 

It is assumed that the junction surface between modules has no deformation. 
Another assumption is that appendages may not be connected to each other; i.e., 
a closed loop may not be formed. However, the user may bypass the Phase I 
computer program and supply the normal modes and modal masses of any Laboratory 
and/or Counterweight to Phase II. 

One feature of the procedure is that the user is permitted to supply con- 
strained substructure modes. These are modes which were obtained for idealizations 
where constraints were employed; for example, in a beam analysis axial extension 
may have been neglected. Constraints may be handled by two different methods. 

In the first, the user supplies modes containing six rectangular coordinates for 
each mass point; however, constrained coordinates are either related by equations 
of constraint or may be zero. By the second method, the user may supply modes 
containing less than six coordinates to describe the motion of each mass point; 
however, he must either identify which if any of the six coordinates are zero or 
he must supply constraint relations which transform the supplied coordinates to 
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Fig. 2.1 General 17-Module Configuration for Modal Synthesis 
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the six rectangular coordinates. As described in Section 2.3.2, only certain 
types of contraint relations may be supplied. This feature was limited by the 
available computer storage space. 

2.2 SELECTION OF SUBSTRUCTURE MODES 

The normal modes of vibration of each of the substructures must first be 
determined fr om conventional lumped-parameter vibration-analysis procedures. 

The required input data includes the natural frequencies, mode shapes, generalized 
masses, and discrete mass matrices associated with each flexible body. As 
indicated in Figure 2.2, if free -free modes of the coupled structure are desired, 
then the free-free modes of Module 1 must be supplied; if the cantilevered modes 
of the coupled structure are desired, the the cantilevered modes of Module 1 must 
be supplied. Module 1 may also be contrained in some directions and not others. 
When obtaining the modes of Module 1, it must be supported in the same way as 
the total structure. Accordingly, the modes of Module 1 may have from 0 to 6 
frequencies which are zero... Cantilevered modes are used for all other modules. 

All of the substructure modes with the exception of the rigid-body modes for 
Module 1 must be orthogonal; therefore the user must ensure that any equal- 
frequency elastic modes have been orthogonalized before use in this program. 

As with any synthesis technique, the procedures discussed here will yield 
an exact representation of the coupled system if all of the substructure modes 
are used. In practice, however, each substructure is approximately represented 
by its lower -frequency modes. Thus the synthesized modes are approximate, with 
the lower-frequency modes generally being most accurate. How well the synthesis 
procedure works is dependent on how many modes are used to represent the sub- 
structures, and how accurately these substructure modes represent the motions of 
the coupled structure. 

The analyst can exercise some control over the accuracy of the synthesis 
through the judicious selection of the substructure modes to be used. Another 
means of improving accuracy is to use so-called "mass-loaded" modes. This means 
that the mass and inertia properties of all influencing substructures are included 
in the prel imi nary idealization used for obtaining the modes of a core module. 

For example, referring to Figure 2.2, the modes of Module 1 would be calculated 
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INNER JUNCTION POINT 
OF MODULE 11 
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Fig. 2.2 Arrangement of Substructures for Modal Synthesis ; 
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with Module 2 through 1J attached as rigid masses, those of 11 with 12 through 17 
attached, those of l4 with 15, 16, and 17 attached, and those of 16 with 17 attached. 
The modes obtained in this manner are closer to those of the coupled system than 
modes computed without mass loading, and thus produce consistently better results 
when used in the synthesis procedure. See Section 3.4 of Volume I for a numerical 
comparison of frequencies using mass-loaded modes and frequencies using conventional 
modes. 

2.3 PREPARATION OF DATA 

2.3.1 Numbering System and Coordinates 

Each module is assigned a module number. Its mass points are numbered 
consecutively beginning with 1 as shown in Figure 2.3. The input modes are 
expressed in terms of these numbers. The program renumbers the mass points 
consecutively for the entire structure, starting with Module 1 and continuing, 
in order, with each higher module number. The coupled structure numbers are 
shown in parenthesis in Figure 2.3. The output for the entire structure is 
expressed in terms of these numbers. 

Each module (except Module l) must attach to a mass point of the previous 

module; however, if there is no mass present at the connection point, an 

artificial mass may be created which has zero mass.* In the example. Module 14 
of Figure 2.3 connects to Mass 3 of Module 11. The attachment point (Point 0) 

of each module is located by its global coordinates (X^, Y^, and Z^). Point 0 

for Module l4 and the global coordinates (X, Y, and Z) are also shown in Figure 2.3. 5 

The global coordinates are also the mean axes (the X axes for the Laboratory and 
the Y axes for the Counterweight) of the Phase II program. 

Since the user's modes are supplied in a local (x, y, z) coordinate system 
■which is not generally parallel to the global coordinate system, the orientation 
of the local coordinates must be specified for each module. This is accomplished 
by the method illustrated in Figure 2.4. As indicated, the local coordinate 
system of a typical module is translated so that its origin will coincide with 
the origin of the global coordinate system. Then points a and b are selected 
as follows : 

* In this case, it is advisable to use the mass-loaded modes discussed in the 
previous subsection. 



MASS POINT 
NUMBER 



Fig. 2.3 Example of 
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Fig. 2.4 Information Required to Specify Orientation of Local Coordinate System 
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• Point a is any point on the positive z axis 

• Point b is any point in the first quadrant of the x, z plane that is 
not on the z axis 

It is usually convenient to select point b as a point on the positive x axis. The 

, Z^) are specified and 
this information is used by the program to orient the local coordinate system. 

2.3*2 Constrained Substructure Modes and Nonrectangular Coordinates 

The program normally requires, as input data, substructure mode shapes 
expressing the motion of each mass point using three rotations and three translations 
expressed in the local coordinate system for that substructure. However, as a 
result of the method used for idealizing the substructure, the user may have modes 
available with less than six degrees of freedom at each mass point. As an example, 
if the substructure is a beam, axial extension may be neglected. In other types 
of constraints some coordinates are dependent on others. In addition, coordinates 
which are not rectangular (for example, polar coordinates) may have been used in 
describing the substructure modes. The program can accommodate modes which do 
not use the usual standard six rectangular coordinates to describe the position 
of each mass provided the N (N ^ 6) coordinates used {q^} at each mass point 1YL 
can be related to the six rectangular coordinates [q^} at fh by the following 
equation: 


global coordinates of points a and b (X , Y , Z and X. , Y, 

& a SI D D 


lf,i ■ ir.flf'l '“> 

The computer must be informed of the value of the N by 6 T matrix Tibi 
at each mass point; the following program features simplify the required input 
data in most cases: 

• AT matrix need only be input once if it applies to an entire module or 
to consecutively numbered mass points within a module. This matrix is 
known as the default T matrix. 

• If all degrees of freedom present at a mass point are expressed in the 
rectangular local coordinate system (e.g., x, y, z for Module 14 in 
Figure 2.3), it is only necessary to indicate which coordinates are absent 
in the mode shapes as a result of zero deformation in that direction. 

The program will construct the T matrices from this information. If all 
coordinates are present, the T matrix will simply be the identity matrix. 
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2.4 PROGRAM LIMITATIONS 

1) Module 1 must be present 

2) Modules 2-11 if present must attach to Module 1 

3) Modules 12 -l4 if present must attach to Module 11 

4) Modules 15, 16 if present must attach to Module 14 

5) Module 17 if present must attach to Module 16 

6) There is a maximum of 30 input-data modes per module 

7) There is a maximum of 100 total input-data modes 

8) There is a maximum of 30 mass points per module 

9) There is a maximum of 100 mass points total 
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2.5 INSTRUCTIONS FOR SETTING UP INPUT DATA DECK FOR PHASE I PROGRAM 

The input deck to the program is organized as follows for each run: a 

TYPE I deck, then several TYPE II decks, followed by a TYPE III deck. The 
user may supply two runs back-to-back. If only one run is supplied the program 
automatically detects the end of the input data deck and terminates execution. 

The TYPE I deck describes which modules are present, the number of modes 
per module, the number of mass points per module, junction points for attached 
modules or appendages, and the gravity constant (G) for the program. A TYPE II 
deck is enclosed for each module and gives the location of the attachment point 
(Xq) and orientation of the local coordinates (X & and X^) for the module, each 
mass point location in local coordinates, the mass matrix in local coordinates 
corresponding to each mass point, the module modes in local coordinates, and 
data identifying the degrees of freedom present at each mass point. The TYPE III 
deck includes the diagonal elements of the modal mass matrix and the substructure 
eigenvalues (frequencies squared). 

Specific instructions are presented below. An example of the input data is 
given in Appendix A-l 

TYPE I Deck 

The TYPE I deck consists of 5 cards. 

Card 1: Seventeen two-column fields for modules 1-17. Place a 1 in the two- 

column field corresponding to a module which is present. 

Card 2: Similar to Card 1. Place the number of modes for each module present 

in the appropriate two-column field. 

Card 3: Similar to Card 1. Place the number of mass points for each module 

present in the appropriate two-column field. 

Card 4 : Sixteen two-column fields for modules 2-17. Place the number of the mass 

point to 'which the module is attached in the appropriate two-column field. 
Use the mass point number of the root module. (E.g., Module 14 would be 
connected to mass point 3 in Figure 2.3.) 

Card 5: Enter the gravity constant G in the first 15 columns. The elements of 

the mass matrix are divided by G. It is read with a FORTRAN format of 
(E 15.8). 
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TYPE II Deck (For Each Module Present) 


Group 1: Enter the global coordinates, (X, Y, Z), of points 0, a and b with 

(X n , Y_, Z.) on the first card, (X , Y , Z ) on the second card and 
(X^, Y fc , Z b ) on the third. Each card has three 15 column fields 
corresponding to the X, Y and Z coordinates respectively. The FORTRAN 
format is (3E 15*8) for each card. 


Group 2: Enter the local (x, y, z) coordinates of mass points 1, 2, 3, etc. for 

the module. Use one mass point per card with FORTRAN format (3E 15-8). 


Group 3: Enter the mass and the 3x3 inertia matrix for each mass point 

1, 2, 3, etc. consecutively. Mass and inertia properties may be 
supplied in mass or weight units, along with the appropriate value of 
"G". The actual mass and inertia properties are input without mass 
loading. Four cards for each mass point must be supplied: 

Card 1: Enter mass in first 15 columns. The FORTRAN format is (E 15.8). 

Card 2: Enter row 1 of inertia matrix. The FORTRAN format is (3E 15.8) 

for each column. 

Card 3: Enter row 2. 

Card 4: Enter row 3. 


Group 4: This group specifies the module mode shapes and the constraint matrices 

(T matrices) discussed in Section 2.3.2. 

f 

(a) Enter a T-matrix card to define the default T-matrix for the 
modules. This matrix will be used for all mass points unless 
overridden. A T-matrix card is a seven-field card. The first 
column contains a 0 if the coordinates used for the mode shapes 
are parallel to the local axes and it is desired that the program 
generate the T matrix. Otherwise, enter a 1 in the first column. 
If a zero is entered in the first column, the next six columns 
contain a 1 if the degree of freedom is present in the mode shape 
and a zero if not. If a 1 is entered in the first column, the 
next six columns should be left blank. 
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(i) (If column 1 is set to l) - Enter N on a card, where N is 
the number of degrees of freedom used to describe the motion 
of the constrained mass point in the mode. Use a FORTRAN 
format of (12). 

(ii) (if column 1 is set to 1) - Enter the NX6 T matrix row by 
row, four numbers to a card, using format (4E 15 . 8 ). 

(b) Enter the coupled-structure mass-point number of the first ma ss 
point at which the default T matrix is to be overriden. Use 50 
if the default T matrix is not to be overridden. Use a FORTRAN 
format of ( 12 ) . 

(c) For each mass point: 

(i) (if override was specified for the current mass point) 

Perform steps (a) and (b) above to define the override T 
matrix for the mass point and the number of the next mass 
point at which the default is to be overridden. Do not enter 
the next mass-point number when overriding the T matrix at 
the last point . 

(ii) Enter module modes in local coordinates. If any rigid -body 
modes are present, these must be the first modes in Module 1. 
The modes are read in as an NXM matrix row by row where N 

is the number of degrees of freedom at the current mass point 
and M is the number of modes for the module. The FORTRAN 
format is (3 (8X, E15.8, IX)). 

Repeat the TYPE II Deck for each module present. 

TYPE III Deck 

Group 1: For each module present enter the diagonal elements of the modal mass 

matrix using format (3 (8X, E15.8, IX)). Begin a new card for each 
module . 

Group 2: Similar to Group 1. Enter the eigenvalues (frequencies squared in rad 2 

per unit time squared). 
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Group 3: Consists of 2 cards. 

Card 1: Enter the number of rigid-body modes in column one. If the 

program is to be used together with the Phase II program, 

6 rigid-body modes should always be used. 

Card 2: Enter the number of eigenvectors to be computed, and the number 

of eigenvalues to be computed. The FORTRAN format is (213) • 

The number of eigenvectors must be less than or equal to the 
number of eigenvalues, and the number of eigenvalues must be 
less than or equal to the total number of substructure modes. 
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2.6 PROGRAM OUTPUT 

The results of this program are the flexible modes of the coupled 
structure in the global coordinate system. The trivial rigid-body modes 
are not printed. The input data for the Phase II program must not contain 
the rigid body modes*; thus the Phase I output modes are in the foirn re- 
quired for the Phase II program. The Phase I output modes are normalized 
to unit modal mass. In addition to the combined- structure mode shapes and 
frequencies, the Phase I program output includes the mass-point locations 
in global coordinates and the mass and inertia matrices in global coordinates 
for each mass point. This information is also required input for the Phase 
II computer program. A sample of Phase I output is presented in Appendix A 2. 

The Phase I and Phase II computer programs may be connected and run 
together by using the proper job control instructions. To accomplish this, 
the output of a Phase I run is placed on a file which is then passed to the 
Phase II run. If two back-to-back Phase I runs are made, two separate files 
are created. The output files of the Phase I program are discussed in 
Appendix A5> and the input files of the Phase II program are discussed in 
Appendix b6. 


* The Phase I input data, however, must contain the rigid-body 
modes of Module I. 
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3.0 PHASE II - DEVELOPMENT OF TIME HISTORY 
3.1 INTRODUCTION 

The Phase II computer program solves the equations of motion developed in 
Volume I by numerical integration and, in this way, develops the time history of 
the system's motion. The entire vehicle under consideration is referred to as 
the Space Station, the main body is called the Laboratory and the other body is 
called the Counterweight. The structure separating these bodies is referred to 
as the Connecting Structure. Both the Laboratory and the Counterweight are 
flexible structures with arbitrary shape and mass distribution. The Connecting 
Structure is also flexible and arbitrary, but its mass is neglected. Figure 3-1 
illustrates the general configuration which may be treated, and Figure 3.2 
illustrates a particular example. 

Problems which may be studied using the simulation include deployment and 
retraction of the Connecting Structure with simultaneous or sequential spin up 
and spin down, the effects of moving rigid masses such as a cargo elevator or 
crew members on board the Laboratory, and the effect of fluid pumped through a 
piping system on the Laboratory. Five control systems are included to study 
various types of controlled maneuvers including mass balancing and wobble damping. 
Various features including control systems and Connecting-Structure characteristics 
have been included in subroutine form so that these items may. be easily replaced 
with minimum disruption to the main program. Special constraint options permit 
the user to rigidize the Laboratory, the Counterweight, or the entire Space 
Station; thus the easier-understood rigid-body results can be compared with the 
more complex flexible-body solutions. 
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Fig. 3.1 Idealization and Coordinate Systems 
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Fig. 3.2 One Example of a Space Station that Can be Studied Using the General Mathematical Model 
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3.2 SYMBOLS 


The analytical symbols used in Section 3 and in Appendix B5 are 
defined below. 


A 

AE 

a 

a 


acc 

acc 

acc 


01C 

^03 MAG 
MAG 


b 


b. 

l 

E,E 


El 

e u)lDB 

e 0DB 

e 0MAX 


err KT ,err KTI 


err LL’ err LLI 

K) 


{Ae i } 


cross sectional area of inside of pipe 

Axial stiffness of beam-type connecting structure 

undeformed distance from Laboratory mass center to 
attachment point of beam-type Connecting Structure 

index number of mass point on Counterweight to 
which Connecting Structure is attached 

magnitude of spin-speed acceleration command 

magnitude of deployment acceleration command 

magnitude of acceleration command 

undeformed distance from Counterweight mass center 
to beam attachment point for beam-type Connecting 
Structure 

inside radius of reservoir on m. 

i 

points at which Laboratory and Counterweight, 
respectively, attach to Connecting Structure 

bending stiffness of beam-type Connecting Structure 

spin speed control error dead band 

attitude angle error dead band 

magnitude of maximum desired attitude angle error 
rate at which torque is applied 

used in both gimbal-angle and gimbal-angle rate 
law (see Figures 3*1^ and 3*15) 

unit vector tangent to pipe centerline along nominal 
direction of fluid velocity ^ ^ 

value of fe.3 at nominal pipe outlet less value at 

1 ( 1 ) 
nominal pipe inlet 


(l) See Footnotes at end of symbol list 
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f JET 

CfiUfi! 

t f * E ) 

f°iJ 

re’’} 

GJ 

fAH.} 

h. 

1 

h. ,h. ,h. 

10 ’ 1 max l min 


I B GIM 


{e } for pipe section on m. which is connected to 
i b ^ 

reservoir on iru . fe^) is a i- wa y s directed outward 

from control volume (see Figure 3*6) 


unit vector along reservoir center line in mass nu 
directed from reservoir base toward reservoir-pipe 
connection ^ (see Figure 3*6) 

jet thrust 


internal force at point on m . , and m“ , expressed 

V ~V ^ 

X and Y axes, respectively 

/V Arf 7 

supplementary forces on m. and m , respectively^ ^ 

1 3. 

force applied by Laboratory on Connecting Structure at 
point E; vector is expressed in :z axes 

force applied by control system on m. or m , 

.. . (1) (2) 1 a 

respectively 

see Equations (3. 4a) and (3.4b) for a through pipe^ 


{G^} for a reservoir pipe 


( 1 ) 


torsional stiffness for beam-type Connecting Structure 


[F(U. )] [e } at nominal pipe outlet less value at 

1 1 (1) (6) 

nominal pipe inlet for pipe on m, 

fluid height in reservoir on m^ 

initial, maximum and minimum values of h^, respectively 

transverse moments of inertia of undeformed Laboratory 
and Counterweight, respectively, used in simplified 
model to determine damping coefficients of beam-type 
Connecting Structure in bending 

moment of inertia of CMG wheel including gimbal about 
transverse wheel axis 

moment of inertia matrices of m. and m , respectively^ ^ 

i a 


ci*], rij] 


moment of inertia matrices of fluid in reservoir pipe 
on through pipe on m respectively 
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t 


i 


i 


s 


J L ,J C 


rjf] 


r aj. ] 

i 


index of mass on Laboratory to which Connecting 
Structure is attached 

index of mass on Laboratory at which sensor and 
CM3 are located 

axial (roll) moments of inertia of undeformed 
Laboratory and Counterweight, respectively used 
in simple model to obtain damping coefficient of 
beam-type Connecting Structure in torsion 

area moment of inertia of the elements dA on the 

n (1) 

fluid boundary about point A i (see Figure 3- 5b) 

l"j\] evaluated for area associated with point B 
on reservoir pipe (see Figure 3.6 a) 

Tj.] at nominal outlet of through pipe less value 

\ . . . . ( 1 ) ( 6 ) 

monimal inlet v ' 


K eXSE 


gain in balance-mass control system multiplying 
filtered sensor acceleration error 


Sese 


Ki 


gain in balance-mass control system multiplying 
balance mass velocity 

gain in wobble control system multiplying error 
in gimbal angle and gimbal-angle rate in gimbal-angle 
and gimbal-angle rate control system respectively 

gain in wobble control system multiplying gambal-angle 
rate 


component 3 of [ }, length of beam-type Connecting 
Structure 


*030 


W,U 0 ) 


update value of i specified for use in position 
command routine when deployment or retraction is 
commanded 

deformed and undeformed vectors, respectively, from 

Connecting-Structure attachment point on Counterweight 

( 5 ) 

to attachment point on Laboratory w 

minimum value of H used in evaluating stiffness and 
damping coefficients for beam-type Connecting Structure 


0 min 
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f 

m. 

1 

P 

m. 

1 


mom 


CMG 


n,n 


n 


CMG 


P,P 




in list of input data, these values are masses of 
Laboratory and Counterweight, respectively, used in 
beam-type Connecting Structure to establish damping 
coefficients 
"t h 

mass of i lumped mass point of Laboratory or Counter- 
weight, respectively 

mass of fluid in pipe on iru 

mass of fluid in reservoir pipe on m 

momentum of CMG used for wobble control 

number of mass points on Laboratory and Counterweight , 
respectively 

value is 1 when CMG is nominally controlling wobble 
(gimbal-angle mode) otherwise value is 0 (CMG gimbal 
angle -rate mode ) 

number of flexible modes used for Laboratory and 
Counterweight, respectively 

linear deflection of cm of m. or m , respectively^ ^ 

x a 


{r} 

vector from origin of Z axes 

to 

origin of X axes 

[r} 

vector from origin of X axes 

to 

origin of Y axes 

r>u 

«„> 

undeformed value of f R } 



{AR} 

{r} - Cr o ) 



{R.} 

vector from origin of Z axes 

cm 

of m. 

l 

fr.} 

vector from origin of X axes 

to 

cm of m. when 

l 

~ ( 3 ) 

Laboratory is undeformed 


l 

^a J 

vector from origin of Y axes 

to 

cm of m when 
a 


Counterweight is undeformed 

■ 


( Sl ) 

see Equation (3. 3b) ^ 



fs. p } 

1 

value of {S.} for reservoir pipe on m 

i l 


(3) 

(4) 


d 
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(s } 




fi] 


t 



t i 03 U 

t qLIM 


t culCU 

ft r } 


(t^, } 


fu },{u.(o)} 

J <3 


fu.} 

1 

f Au. } 

1 

V 

v u 

Vel Jl03 MAX 

Vel MAX 

XJ1 


vector from cm of m to Connecting -Structure attachment 


point on Laboratory 


- ( 1 ) 


vector from cm of m to Connecting-Structure attachment 

~ ( 2 ) 

point on Counterweight 
time 

time of update 

time at which fluid velocity is updated 


time at which £ is updated 
limiting torque of CMG gimbal motor 


time at which angular- velocity command is updated 

vector from cm of m t to point on m , where internal 

(\v J 1 

load is required ' ' 

vector from cm of in", to point on m~ , where internal 
load is required 

position update of moving mass 


position and initial position, respectively, of 

(3) 

moving mass when Laboratory is undeformed 

(j = 1 designates elevator, j = 2 designates balance 

mass ) 


third component of fu } 

3 

vector from cm of m. to center line of pipe on 


m. 


(1) 


fu } at nominal pipe outlet lass value at nominal 
pipe inlet ^ ^ 


fu^} at points R and B respectively (see Figure 3.6) 
fluid velocity in pipe 


command update of fluid velocity in pipe 

maximum value of 

maximum commanded velocity 

jet moment arm (see Figure 3*8) 
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mean axes for Laboratory- 

body axes in iru with origin at cm of m^ 


coordinates of accelerometer in X s axes 


^ X SEWS/CM^’^ X SENS/CM^ commanded vector and actual vector, respectively, 

from Space Station center of mass to sensor (accelerometer) 

( 3 ) 

for use in balance mass control system 
Y mean axes for Counterweight 




body axes in m with origin at cm of m 
8 . 8 

see Section 3-3.^ 
inertial coordinate system 


Connecting-Structure coordinate system fixed in m 


[e.},{0 } 

i a 


CMG gimbal angle 

Euler angles orienting the £ axes with respect to 
the Y~ axes 

cross-product function 

ratio of actual to critical damping coefficient of 
beam-type Connecting structure for axial, bending, 
and torsional vibration modes , respectively 

Euler angles locating X axes with respect to Z axes 

deformation of Connecting Structure, 

Euler angles orienting Y axes with respect to X axes 

fn) - jo 

angular deformation of m. and m , respectively 

i 8 

_ .th 

mass of j moving mass 


^1EST 


{§],{§} 


estimated mass of elevator for use in balance mass 
control system 

model coordinates for Laboratory and Counterweight, 
respectively 
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[n ( )] 


coordinate transformation matrix 


P 

t bal 


f T f 1 f T“ 

1 V j ’ 1 V 


} 


{$*} 


CD , U) 

1C’ 1CU 


(i) 

1C ML 




} 


mass density of fluid 

filter time constant used in balance mass control 
system (set to zero to bypass filter) 

torque exerted by Laboratory on Connecting Structure 
at point E; vector is expressed in z axes 

internal torque at point on m. t and m- , , expressed 

JL' l' ^ ^ 

in X and Y' axes, respectively 

Euler angles expressing angular deformation of mi- 
relative to m~ in z axes 

IN/ 

commanded angular velocity and angular velocity 
update, respectively 

minimum value of U) at which balance mass control 

XU 

system is in operation 

angular velocity of X and Y axes, resepectively ^ ^ 
angular velocity of X^ axes ^ 


Footnotes : 

(1) The matrix form is expressed in the X 1 axes for the symbol pertaining 
to the Laboratory 

(2) The matrix form is expressed in the Y a axes for the symbol pertaining 
to the Counterweight 

( 3 ) The matrix form is expressed in the X axes for the symbol pertaining 
to the Laboratory 

(4) The matrix form is expressed in the Y axes for the symbol pertaining 
to the Counterweight 

( 5 ) The components of the vector are expressed in the z axes 

(6) Nominal pipe outlet and inlet refers to outlet and inlet, respectively, 
when fluid velocity V is positive 
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3.3 COORDINATES USED IN THE PROBLEM 
3.3.1 Coordinate Systems 

Figure 3.1 illustrates the coordinates used in the problem. The axes of 
each coordinate system used are always denoted as the 1, 2, and 3 axes; for 
example, the Z coordinate -system axes are Z^, Z^, and Z^. The main coordinate 
systems are the Z axes, the X axes and the Y axes. The Z axes are an inertially 
fixed reference coordinate system. The X axes are a system of mean axes moving 
with the average motion of the Laboratory and are used as the reference system 
for elastic motion of the Laboratory. The Y axes are similarly a set of mean 
axes for the Counterweight. The X and _Y coordinate systems are also the global 
coordinate systems in which the modes of the Laboratory and Counterweight, 
respectively, have been computed. 

The axes fixed in each mass mu of the Laboratory are the X 1 axes with their 
origin at the mass center of nr. The Y 2, axes are similarly a set of axes 
fixed in a typical Counterweight mass m . 

cl 

The axes must be selected so that they meet the following orientation 
specifications . when the system is undeformed . First, the X axes are selected as 
any convenient reference system for the Laboratory. AL1 of the X axes must be 
parallel to X. The Y coordinate system is rotated l80° about the number 1 axis 
with respect to the X coordinate system. Thus, the X and Y coordinates are 
parallel; however the Y g and Y^ axes are opposite in direction to the X g and X^ 
axes, respectively. The Y 2, axes are parallel to the Y axes. These parallel 
requirements on the nominal positions of the coordinate systems do not limit the 
generality of the program since it is always possible to select the coordinates 
with the orientations required. In fact, if mass distribution and modal ^ 
properties of a structure are available in other coordinate systems, the Phase I 
program may be used to transform these properties into the required coordinate 
systems. 


3.3.2 Euler Angles 

Each set of Euler angles is designated by a 1, 2, and a 3 angle; for 
example, the set of angles used for the coordinate transformation from the Z 
to X axes of Figure 3.1 are called [y] and the specific angles are y-^ y 2 ’ 

Y . The order of rotation is always as follows: the number 1 angle about 

the 1 axis, the number 2 angle about the carried number 2 axis, and the number 
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3 angle about the twice-carried number 3 axis. This is illustrated in Figure 
3-3 for the angles {v}. In the figure the X axes have been translated so that 
their origin coincides with that of the Z axes for the purpose of illustrating 
the method of orienting the X axes with respect to the Z axes. For brevity 
{y} is simply described as the "set of Euler angles orienting the X axes 
relative to the Z axes." Similar brief descriptions will be used to define the 
other Euler-angle sets. 


3-3-3 Main Coordinates Used in the Problem 


As indicated in Figure 3.1, { R} and the Euler Angles f y } locate the X 
axes in space, and {R} and {T|} locate the Y axes with respect to the X axes. 

The elastic linear and angular deflections are {q^} and {0.}, respectively, for 
a typical mass point rru (i = 1, ..., n) on the Laboratory and {q^} and { 0 & } for 
a typical mass point in (a = 1, . . . , n) on the Counterweight. These elastic 

cl 

coordinates were linearized whenever there was an analytical or computational 
advantage to do so; however, the rigid-body coordinates {R}, {y}, {R}, and {Tl} 
are not linearized so that large rigid-body motions may be studied. Instead of 
Euler-angle rates ({y} and { T) } ) , angular velocities ({cu X } and {ou Y }) of the X 
and Y axes, respectively, are used as coordinates. To improve numerical accuracy, 
the following perturbation parameters were used as integration coordinates in 
place of {R} and {Tl}: 


{AR} = {R} - [R o ] 


(3.1) 


f'Q*} = [T]} - 0 

(o 


(3.2) 


The derivatives of {r} , {y} , {u) X } , the {qj’s, the {6.}*s, [AR] , {!}*}, 

{m }, the { } ’ s and the { © a } 1 s are obtained in the program and numerically 
integrated. In addition certain control-system variables are also integrated. 
The integrated variables are given in Equation (3.5). 


3-3.4 Coordinates Used to Describe the Connecting-Structure Deformations 
It is assumed that the Connecting-Structure is fastened to m on the 

cl 

Counterweight and m. on the Laboratory. (See Figure 3-1- ) The z axes which 
are fixed to m & at point E are a convenient set of coordinates for expressing 
the Connecting-Structure deformations. The origin of these axes, point E, is 
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/ 

/ 

/ 

( 2 ) 

SECOND ROTATION 



THIRD ROTATION 


Fig. 3.3 Typical Set of Euler Angles 
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located with respect to the center of mass of m by specifying the vector {s}. 

cl 

rv/ 

The orientation of the £ axes is specified in the same manner as discussed in 
Section 2.3.1; i.e., by specifying the vectors [Y^. and {yS. } in the Y^- 

^a j J i^b) ~ 

coordinate system which is parallel to the Y^ coordinate system but has its 

<■**»* 

origin at E. These vectors locate point (a) which is any point on the positive 
axis and point (b) which is any point, in the first quadrant of the z.^ plane. 

Point E is a reference point on nu on the Laboratory. The vector {l}, from 

E to E, may be thought of as the length vector of the Connecting-Structure. [£} 
is expressed in the £ axes and is composed of the undeformed length {£ Q } plus the 
deformation {6}. Luring deployment and retraction maneuvers, [j2 } is a specified 
function of time. When there is no deformation, a coordinate system z may be 
imagined as fixed in rru and parallel to z. During deformation, the Euler angles 

expressing the angular deformation of the z axes relative to the z axes are 
Thus, the Connecting-Structure deformation is specified by the vectors {6} and 
{ 0 *}. 
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3.4 ENGINEERING CONSIDERATIONS WHEN PREPARING DATA 

The input data deck is to be constructed in accordance with the procedure 
described in Section 3*5. Special engineering considerations are discussed in 
this section. 

3.4.1 Constraint Options 

Constraint options may be used to run the program with portions of the 
vehicle rigid even though the available data deck was prepared for a flexible 
vehicle. This enables the analyst to compare the more complex flexible-vehicle 
results with simpler rigid-body results. It is possible to: 

• rigidize the Laboratory 

• rigidize the Counterweight 

• rigidize the entire Space Station 

• rigidize both the Laboratory and the Counterweight while maintaining 
a flexible Connecting Structure. 

3.4.2 Basic Input Data for the Laboratory 

The following basic input data is required to define the mass and flexi- 
bility characteristics of the laboratory: 

• Coordinates of the center of mass A. of each mass point m. in the X 

i i ~ 

coordinate when there is no deformation. (See Figure 3-l). 

• Mass irn of each mass point and the inertia matrix [1^] in local 
coordinates X 

• Elastic frequencies and mode shapes relative to the X coordinate system. 

The free-free modes are required and the six rigid-body modes must be 
omitted. The modes must be orthogonal. Each mode must contain the three 
translational coordinates and the three rotational coordinates, in sequence, 
for each mass point. If the moment of inertia, at a mass point is neglected 
and there is no torque at that mass point, the rotational coordinates are 
not used; thus, the' corresponding field on the data card may be left blank 
if desired. 

As indicated in Section 2.6, it is possible to generate all of the above 
information with the Phase I computer program, and, if desired, connect the 
two programs so that the information is passed directly to the Phase II pro- 
gram. This information may also be generated elsewhere and input directly 
into the Phase II computer program. 
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3.4.3 Basic Input Data for the Counterweight 

The basic input data for the Counterweight is analogous to that described 1 
above for the Laboratory, and the above comments are also applicable to the 
Counterweight data. The modes and inertia matrices for the Counterweight must 
be expressed in the Y and Y 3 axes, respectively (see Figure 3.l). 

3.4.4 Initial Conditions and Mean Axes 

While the positions of the X axes for the Laboratory are arbitrary when 
there is no deformation, when the structure is deformed the coordinates move 
with the mean, or rigid-body, motion of the Laboratory. The program computes 
the position and rate of movement of these axes and prints out this information 
as {R}, { Y }, {R} , [w }. When the system is initially deformed it would be a 
burden for the user to average the deformation and deformation rates and compute 
the initial position and rate of movement of the mean axes. To avoid this bur- 
den, the user may supply the initial deformation data ({cu}, { 0 i ) , {q.}, and 

[0. }) in terms of any neighboring axes (defined by the initial values of {R}, 

± 

[y], {R} and [oj }) and the program will compute the initial position and rate of 
the mean axes and convert the initial deformation data to these axes. For this 
reason the values of {R} , { Y } , {R} , {</}, {q^, {0J, {q^, and {0J (for i=l, 

..., n) that are printed out at t=0 may not agree with the values that are input 
to the program. 

The Counterweight is located relative to the Laboratory, at t=0, by speci- 
fying the initial undeformed location of the Connecting Structure its deformation, 

({6} and {0*}) and the deformation rates {6} and {0*}. The elastic coordinates 

• • 

of the Counterweight {q^}, {0J and their rates {q^, and {0 } for a=l, ..., n 
are initially described relative to the Y axes. The Y axes are a set of mean 
axes for the Counterweight, and, as is the case for the X axes, the program will 
initially recompute the location of the Y axes and transform the elastic coordi- 
nates. 

3-4.5 Damping in the Laboratory and Counterweight Struc tures 

Modal damping is assumed. The ratio of the damping coefficient to the 
critical damping coefficient is specified for each mode. 

3.4.6 Fluid Motion 

The dynamic effect of fluid being pumped through a piping system may be 
studied using this program. It is assumed that the pipes terminate in two 
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reservoirs, one called the emptying reservoir and the other called the filling 
reservoir . (See Figure 3.4 a. ) These reservoirs are emptying and filling, 
respectively, when the fluid velocity V is positive; however if the fluid 
reverses direction, the emptying reservoir would be filling and vise versa so 
that the names of the reservoirs are reference designations only. 

The first step in preparing the fluids data is to separate the fluid 
system into its elements. To do this, control surfaces are drawn about each 
structural mass point associated with fluid motion as shown in Figure 3-4A. 

The fluid contained within each control surface is assumed to vibrate with the 
associated mass point as shown in Figure 3-4B, in exaggerated form; thus, the 
selection of the control surface boundaries involves judgement of the user. 

Each control surface will contain either a through pipe (as shown in Figure 3-5) 
or a reservoir and reservoir pipe (as shown in Figure 3-6). 

3. 4. 6.1 Through Pipe 

For each through pipe, the user must specify the mass distribution and geo- 
metrical properties of the fluid in the through pipe. Some of these terms, such 
as S^, require discussion. is defined as follows: 



(3.3a) 


where dm is a differential fluid mass on nn , u* is the vector from A^ to dm and 

the integral is taken over all of the fluid in the through pipe. Alternately, 

— > 

the following formula may be used to evaluate S^: 



(3.3b) 

\ 


where ds is the differential arc length along the centerline as shown in Figure 
3.5, u is the value of u' on the centerline, and the integral is taken from the 
nominal inlet to the nominal cutlet. It may also be helpful to notice that 
is simply equal to the mass of the fluid times the vector from A^ to the cm of 
the fluid. 

Another vector required as input data is : 

' OVrLfT 

- J 2? x t/u 

AfV+H/V* i 
t/VCST 


(3.4a) 
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Note: Only two intermediate masses containing 
through pipes are shown; however, the 
analysis was done for an arbitrary number 
of masses 



RESERVOIR 

NOMINALLY 

FILLING 


(A) Undeformed System 


CONTROL VOLUMES ARE FIXED 
TO STRUCTURAL MASSES 



(B) Deformed System 


Fig. 3.4 Fluid Control -Volume Idealization Used for Computing 
Rates of Change of Linear and Angular Momenta 







INLET (FOR V > 0) 


(A) FLUID CONTAINED WITHIN CONTROL SURFACE 


INNER PIPE SURFACE 



(B) DETAIL OF FLUID CONTROL-SUR FACE INTERSECTION 


NOMINAL 
OUTLET 
(FOR V > 0) 


Fig. 3.5 Mass Containing a Through Pipe 
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G^ can be evaluated by expanding u into components or by noting that G^ is the 

following function of the area swept out by u (see Figure 3»7) : 

I/O Af /A/4 t- 

ourcer 

(7 =2-/ ^ ^ A l> (3.4b) 

A/0 M /A/4 t. 

ia/l e r 

where e^ is the unit normal vector which is perpendicular to the plane of u and 
du and directed along u x du. When the surface swept out by u is a plane, is 
constant and the evaluation of Equation (3.4b) is particularly simple. 

The matrix [Jl] is defined as the moment-of-inertia matrix of the elements 
dA^ shown in Figure 3.5b about point A^ in X^ axes. As an approximation, the 

moment of inertia of the projected area shown in the illustration could be used. 

The diagonal elements of [J^] are the moments of inertia and the off-diagonal 
elements are the negatives of the cross products of inertia. 

3. 4. 6.2 Reservoirs and Reservoir-Pipes 

The reservoir is assumed to be cylindrical with radius b^. The pipe-reservoir 
connection is assumed to be along the reservoir axis. The user specifies the 
maximum and minimum fluid heights in each reservoir. If either of these limits 
is reached, the pump is assumed to stop, and no further fluid flow occurs in the 
problem. Input data for the pump command subroutine, which specifies the fluid 
velocity V as a function of time, is discussed in Section 3-4. 9-2. 2. 

Much of the input data for each of the two masses containing a reservoir 
and reservoir pipe is similar to that described in the previous section for a 

P ~ *P f ~ * 

through pipe. For example, [1^] and S^ denote the values of [i^j and S^ for the 

reservoir pipe shown in Figure 3-6A. 

3.4.7 Moving Masses 

There may be up to eight moving masses present on the Laboratory. These are 
designated p^, p 2 , ..., p,g. The subroutines that specify the motion of these 
masses are designed so that they may easily be replaced by the user. A subroutine, 
discussed in Section 3-4. 9.2.1, is supplied with the program so that the user may 
specify a variety of three-dimensional motions for each mass other than p^. p,^ 
is a balance mass which is designed to balance the motion of an elevator along 
axis X 3< The motion of p^ is determined by the balance mass control system 
which is discussed in Section 3.4. 9*2. 6. 

1 ) 
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The motion due to structural deformation is added to the rigid-body motion 
of each moving mass. In order to determine this structural deformation when 
the moving mass is between structural mass points (nodes), an averaging pro- 
cedure is used in the program. The elastic motion of the moving mass is 
approximated as the weighted average of the motion of those nodes that (in the 
user's judgement) influence the motion of the moving mass. The weighted 
average has the following characteristics: 

• As a moving mass approaches an influencing node, its motion will 
approach the motion of that node . 

• The greater the distance a moving mass is from a node, the smaller 
the influence of that node in the averaging procedure. 

The set of influencing nodes could differ for each mass point and could 
vary with time as the moving mass moves through different regions of the Labora- 
tory. This set of node numbers should be programmed into Subroutine LOOK 
(described in Appendix B5) in table look-up form; i.e., as tables of node numbers 
versus time. In the version of LOOK supplied with the program, the influencing 
nodes are 1 through 8 for each moving mass point and remain constant with time. 

3.4.8 Internal Loads 

The facility is provided to compute internal structural loads between 
modules or at any other point where the structure (Laboratory or Counterweight) 
can be separated into two free bodies. This is accomplished by separating the 
structure at the load-computation point and considering one portion of the 
structure as a free body. The load applied by the free body is then computed by 
taking the resultant of all of the inertia and applied loads on the free body. 
Since Connecting-Structure loads are included in the equations as external loads, 
the remaining portion of the Laboratory (or Counterweight ) may also be used as 
the free body. Of course, results which are opposite in sign would then be 
obtained. The index numbers of all mass points in the free body must be supplied 
as input data in addition to the coordinate vector {t^,} ( {t-, } on the Counter- 
weight) to the load-computation point, {t , ) is stated in the local X^ coordi- 
nate system which is fixed in m ,, the mass containing the load-computation point. 

Xj 

The forces {f* ,} and torques {V for the Laboratory are expressed in the X 
axes, and the forces [f-,] and torques {t , } for "the Counterweight are expressed 
in the Y axes. If the option to compute structural loads is exercised, the 

r^/ 

program will also print the structural loads applied by the Connecting Structure 
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to the Laboratory (projected onto X axes) and the Counterweight (projected onto 
Y axes). 

The internal structural loads are neither computed nor printed at t=0 for 
the reason discussed in Section 4.5 of Volume I. 

3.4.9 Data for Replaceable Subroutines 

The physical properties of the Connecting Structure, the control systems, 
the motion of the moving masses, and the externally applied forces have been 
programmed in subroutines which were designed to be easily replaced so that 
different classes of problems could be studied. Details on replacing these 
subroutines are provided in Appendix B5. The input data required for these 
subroutines will now be described. 

3. 4. 9-1 Beam- Type Connecting Structure 

The Connecting-Structure subroutine computes the force and torque applied 
by the Laboratory on the Connecting Structure in the coordinate system. 

The Connecting Structure is assumed massless, and the program uses equilibrium 
to compute the loads on the Counterweight. The subroutine supplied with the 
program contains the equations for a circular beam fixed at each end. During 
deployment or retraction maneuvers it is simply assumed that the undefomed beam 
length increases or decreases, and the beam stiffness and damping coefficients 
are recomputed every cycle based on the updated beam length. Some of the stiff- 
ness coefficients go to infinity as the beam length approaches zero. While 
this occurs in theory, in practice fixed supports are never perfectly rigid; 
thus there is always some flexibility. To approximately account for this 
flexibility and to avoid the numerical problems associated with high stiffnesses, 
a minimum length, min , is input to the program. Whenever i Q is less than 
^0 min’ program uses ^ to compute the stiffness and damping coefficients. 

Damping is based on a simple model of the Space Station. The user inputs 
the desired proportion of critical damping, and the mass properties of a rigid 
Laboratory and Counterweight. The program computes the critical damping of 
this model for each mode, multiplies the critical modal damping coefficients by 
the damping ratios and transforms the modal damping matrix to a physical damping 
matrix. This latter matrix is used in the program. Since the beam is circ ula r, 
the damping coefficients should be the same for in-plane and out-of-plane bending; 
therefore the two transverse moments of inertia of each end body in the simple 
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model must be identical. It is suggested that the average of the two values be 
used for each body. 

3. 4. 9 . 2 General Data for Control Systems 

The control system subroutines compute the forces and torques applied by 
the control system actuators to mass points of the Laboratory and/or the Counter- 
weight. In addition, they may determine the motion of moving rigid masses, the 
undeformed geometry of the Connecting Structure, and the velociy of the fluid 
being pumped through the Laboratory. 

In the control -system subroutines supplied with the program, there are 
sixteen jets. These are arranged to control spin rate and nonspinning vehicle 
attitude for the demonstration Space-Station model discussed in Section 6 of 
Volume I. As indicated in Figure 3*8, there are two clusters of four jets each 
located on mass point 1 (m^ ) of the Laboratory and two clusters of four jets 
located on mass point 5 of the Laboratory. The jet thrust f^^ (which has the 
same value for each jet) and the jet moment arm XJ1 are input data. The jet 
assignments are indicated in the table on Figure 3-8. Twice as many jets are 
used for spin control about axis 1 than for attitude control. 

A sensor package may be located at any mass point i on the Laboratory. The 

s 

sensor is assumed to measure the angular position, angular velocity, and angular 

acceleration as well as the linear acceleration of mass i . 

s 

The controls subroutine group consists of command routines and control 
routines. A command routine is used to specify a certain motion, whereas a 
controls routine is used to achieve that motion. For example, if it is desired 
to vary the vehicle’s spin rate at a prespecified function of time, a command 
routine would be used to generate the desired spin speed as a function of time. 

A controls routine might then fire jets in order to achieve the commanded angular 
velocity. First, the command routines and then the control routines will be 
described. 

3. 4. 9. 2.1 Position Command (For Deployment and Moving Masses) 

The same law is used to command the undeformed position of the Counterweight 
relative to the Laboratory (i.e., third component of {Xq}), and the position of 
moving masses 1, and 3 through 8 ({IL} with i=l, 3? 4, ..., 8). 

The position-command law is illustrated in Figures 3-9 and 3-10. The input 
data for developing the command consists of the initial position, position update, 
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1 



Fig. 3.8 Jet Configuration on Laboratory Core Module 
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ACCELERATION 



Figure 3.9 Position Command Law for Case Where Velocity Command Archieves 
Maximum Velocity 


3-2! 



Fig. 3.10 Position Command Law for Case Where Velocity Command Does Not 
Achieve Maximum Velocity 
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update time, magnitude of acceleration, and maximum velocity. For the case of 
moving masses, a different set of input data may be used for each axis. In 
Figure 3*9 "the motion begins at the update time when a constant acceleration 
occurs until the maximum velocity is achieved. Then cruising occurs at the 
maximum velocity until deceleration occurs at the negative of the acceleration 
magnitude . 


If the position update is close to the initial position, the maximum 
velocity is never achieved. In this case the command law is illustrated in 
Figure 3*10* 


In the event that the position update is less than the initial position 
the negatives of the motions illustrated in Figures 3.9 and 3.10 occur. 


The user must always select a numerical integration time interval that is 
less than the time increment t^ illustrated in Figures 3*9 and 3 . 10 . *01 may 
be computed as follows : 



smaller of 


maximum velocity 
acceleration magnitude 

V | position update - initial position | 
acceleration magnitude 


The acceleration magntude must not be set to zero. If it is desired to 
maintain a constant position command throughout the run, the user should input 
an update time which is greater than the final time designating the end of the 
run. 

3. 4.9.2. 2 Velocity Command (For Spin Speed and Fluid Velocity) 

I 

The same law is used to command the spin speed of the Space Station and 
the fluid velocity. This law is illustrated in Figure 3*11. The required input 
data consists of the initial velocity, the velocity update, the update time, and 
the acceleration magnitude. 


The selected numerical-integration time interval must always be less than 
t Q i of Figure 3.11. t Q1 may be computed as follows: 

' £ _ |update velocity - initial velocity| 

01 acceleration command 


Another precaution that must be taken is never to set the acceleration 
command to zero. If it is. desired to hold the velocity command constant, the 
user should set the update time greater than the final (end-of-run) time. 
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Fig. 3.11 Velocity Command Law 
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3. 4. 9. 2. 3 Attitude Control 

The attitude control system is designed to drive the attitude angles {y} 
to zero for the case of the nonrotating Space Station. The same control law 
is used for each axis. This law is illustrated in Figure 3-12. 

3. 4.9.2. 4 Spin Control 

A hysteresis type spin control law is used. This law is illustrated 
in Figure 3«13* 

3. 4. 9. 2. 5 Wobble Control 

A CMG control law known as the 90° h-lag law is used to eliminate gyro- 
scopic wobble from the rotating Space Station.* The CMG is assumed to be located 

on the same mass (m. ) as the sensor package. There are two modes of control, 

s 

gimbal- angle control and gimbal- angle-rate control. The CMG angle law is used 
for wobble control. When the wobble drops below a predetermined threshold, 
wobble control is ended. The CMG angle-rate law is then used to rotate the 
gimbal relative to the Space Station at the negative of component 1 of the Space- 
Station angular velocity; thus the gimbal angular velocity is approximately zero 
in inertial space. The CMG gimbal angle law and the CMG angle-rate law are 
illustrated in Figures 3-l4 and 3-15, respectively. 

3. 4. 9 . 2 . 6 Balance Mass Control 

The balance mass control system is used to balance mass shifts along the 
long axis of the Space Station (i.e., the axis of Figure 3«l). For example, 
mass shifts would occur due to the motion of a freight elevator. An accelero- 
meter senses the acceleration at a point on mass i . The control system moves 

s 

the balance mass until the centrifugal acceleration at the sensor is equal to 
a commanded value. The commanded vector from the Space-Station center of mass 

to the accelerometer is input data. The user should include the 

elastic deformation of the structure when determining this vector. One method 
of determining this deformation is discussed in Section 3.4.10. A block 
diagram of the control system is illustrated in Figure 3.l6. 

3.4.9. 3 Supplementary Loads 

The user may program any set of loads as a function of time on the Labora- 
tory and the Counterweight into subroutine SUPPLM. The equation for the result- 


* The spin speed must be greater than .02 rad/sec for the wobble control 

I 

system to turn on. 
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Fig. 3.13 Velocity Control 
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Fig. 3.14 CMG Gimbal Angie Control System 
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Fig. 3.15 CMG Gimbal Angie Rate Control System 


£ 
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Fig. 3.16 Space-Station Center-of-Mas$ Position Control System for Flexible 
Vehicle Idealization 
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ant force {f 1 } and the resultant torque {t 1 } on each mass point m. of the Labora- 
s S ^ ^ 1 

tory must he programmed. Also, the resultant force {r } and {t } on each mass 
__ s s 

point m of the Counterweight must he programmed. 

cl 

3.4.10 Determining Quiescent Deflections 

When the Space Station is in a state of pure spin with no vibration, con- 
stant elastic deflections will occur because of the centrifugal forces. This 
state of motion is known as the quiescent state . In most rotating problems the 
user will want to set the initial deflections equal to the quiescent deflections 
or equal to some perturbation about these values. One way of determining the 
quiescent deflections is to make a short preliminary run with of critical 
damping in the Laboratory, the Counterweight, and the Connecting Structure. The 
deflections will quickly damp out to their quiescent values. 
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3.5 PROGRAM LIMITATIONS 

1) A maximum of 20 modes may be used to represent the elastic motion 
of the Laboratory 

2 ) A maximum of 20 modes may be used to represent the elastic motion 
of the Counterweight 

3) There is a maximum of 100 mass points on the Laboratory 

4) There is a maximum of 100 mass points on the Counterweight 

5) There is a maximum of 8 moving masses present 

6) A maximum of 30 masses may contain pipe segments 

7) There is a maximum of l£ points on the Laboratory at which internal loads 
may be printed 

8) There is a maximum of l6 points on the Counterweight at which internal 
loads may be printed 
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3.6 INSTRUCTIONS FOR SETTING UP INPUT-DATA DECK FOR PHASE II PROGRAM 


TYPE I 


Initial Card. - To be supplied at the beginning of the group of multiple 
runs. 

Supply the number of runs to be made, plotting clue [(212) format] 
plotting clue = /0 if no plotting for any run 


{ 1 
ll i 


if at least 1 run has plotting 
TYPE II - To be supplied for each run 

PART 1 - General data for the inn 
Group a 

Card 1 - final time, plot scale [(2E15.8) format] 

(program disregards plot scale if there is no plotting) 
plot scale = number of seconds/inch along abscissa 
Card 2 - plot clue for the run, fluid clue, number or pipes, 

number of moving masses [ (413) format] 
plot clue = /0 no plotting 
' 1 plotting 

fluid clue = i 0 no fluids present 
’ 1 fluids present 

Card 3 - (insert blank card if no moving masses are present) 

enter clues for each moving mass ( 1 - 8 ) [( 813 ) format] 
each clue = (0 moving mass not present 
(l moving mass present 

Card 4 - (insert blank card if no fluids are present) 

mass point index numbers of the emptying and filling 
reservoirs, respectively [( 213 ) format] 

Card 5 - (Card 6 if needed) (insert 1 blank card if no pipes are 

present) mass point index number of each pipe present 
[(1513) format] 

Second card needed if more than 15 pipes are present 

Group b 

(Skip to group C if there is no plotting on this run) 

Card 1 - enter clues for the 13 plotting groups [( 13 II) format] 
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Card 2 


Group c 

Card 1 


KLUE = 


enter 1 if the group is to be plotted, 0 if not. Choose 
7 groups maximum. If plotting is selected on a run {u) X } 
is always plotted first, followed by the plotting groups 
selected on this card. 

The 13 optional groups are: 


Group No. Variables Plotted 


1 

(y3 

2 

{6} 

3 

{**} 

4 

{q 3 for some mass no. a 

cl 

5 

[6^3 for some mass no. b 

6 

[q c 3 for same mass no. c 

7 

{6^} for some mass no. d 

8 

{q g 3 for some mass no. e 

9 

£ o3’ £ o3’ £ o3 

10 

{u 1 } 3 - {u 1 (o)} 3 , [u 2 } 3 - {u 2 (o)}. 


^• X SENS/ CM^ 3 " ^SENS/CM^ 

11 

• 

V, V, hu for emptying reservoir 

12 

[f i for some mass no. r 

13 

f 

{ t 3 for some mass no. s 

Each group consists 

of three plots, one above the other. 


when vectors are plotted, the first component is always 
plotted on top, followed by the second component, with 
the third component on the bottom. 

- enter the mass -point index numbers a, b, c, d, e, r, s 
as above [( 713 ) foimat] 

(insert blank card if none are selected for plotting) 

- enter clues KLUE, IFLAB, IFCW, IFCN [(4l4) format] 

KLUE denotes the type of problem being run 

I I flexible Space Station 
2 rigid Laboratory 

3 rigid Counterweight 

4 rigid Space Station 

5 rigid Laboratory and Counterweight 



3-41 


IFIAB enables the user to run a rigid-Laboratory problem using a 
(slightly modified) flexible -Laboratory deck. 


1 is always used when the Laboratory is flexible; also 
when the Laboratory is rigid and input data for the 
flexible Laboratory is present in the deck, (in this 
case the data is read in but disregarded). 

0 is used when elastic displacements and velocities of 
the Laboratory({q i ], £ 3 , { 3 3 {6^}) are not present 
in the deck. Also, when IFLAB = 0 and the Phase I 
input tape is not present, the Laboratory modal data 
(n umb er of modes, modal masses, frequencies, and mode 
shapes) must be absent from the deck. 

IFCW - same as IFLAB except that it applies to the Counterweight 
input data. 

IFCN - if IFCN = 0, no Connecting Structure data is read in; use 

IFCN=0 only when the Space Station is rigid (see description 
of Part 8). Otherwise set IFCN = 1. 

Card 2 - enter IFDEPL, IFCON, IFSUPP [(314) foimat] 

IFDEPL = ( 0 no deployment 
j 1 deployment 

IFCON = j 0 no controls 

1 1 controls present 

IFSUPP = ( 0 no supplementary forces 
) 1 supp. forces present 



Type 1 controls input (skip this portion if no controls 
are present) 


Card _ l - i g j f JET 


(13, E15.8) format 

Card 2 - Counterweight command input 

X 03U’ t X03U» aCC 403MAG’ Vel 403MAX 
(4E15.8) format 



Card 3 - 


Card 4 - 


Card 5 - 


Card 6 - 


Card 7 - 
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spin command input 

w (initial value), t , nTT , acc , 
1C '* (DlCU* ujic 

(3E15.8) format 

spin control input 

e a)IDB 

(E15.8) format 
attitude control input 

e 8DB’ e 0MAX’ XJ1 
(3E15.8) format 
wobble control input 
P g j 1^, 

(4E15.8) format 

wobble control input (continued) 

*13 > m ° m CMG’ LgIM 


( 3E15 • 8 ) format 
PART 3 - Internal loads 

(if no internal load computations are required, use blank card and 
skip to Part 4) 

Group a - (one card) 

Card 1 - enter clues KLUCS, KHJLAB, KLUCW, N, NB 
( 512 ) format 

use KLUCS =1 if structural loads applied by the Connecting 
Structure are computed. Use 0 otherwise 
use KLULAB =1 if internal loads on the Laboratory are 
to be computed. Use 0 otherwise 
use KLUCW =1 if internal loads on the Counterweight 
are to be computed. Use 0 otherwise 
N = number of load- computation points on Laboratory 

NB = number of load- computation points on 

Counterweight 

Group b - input for each load computation point on Laboratory (skip to 
Group c if N = 0) 

f 

Card 1 - i (the load computation mass-point index number), n umb er of 

mass points comprising the free body [(212) format] 

Card 2 - the load coordinates [t^, } 

(3E15.8) format 
remaining cards (10 maximum) 

the index numbers of the mass points comprising the free 
body, 10 per card 
(1012) format 
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Group c - input for each load computation point on the Counterweight 
(skip to part 4 if NB=0) 
construct similarly to Group b 
PART 4 - Laboratory input data 
Group a - (one card) 

Card 1 - IPHS1, JUMP 

(214) format 

IPHS1 = (0 no Phase I tape provided 

11 Phase I tape provided 

JUMP = 1 is used to set the elastic displacements and 
velocities = 0 without reading them in. 

Group b 

(If IPHS1 = 1, skip to group d) 

Subgroup 1 - (one card) 

Card 1 - enter the number of mass points on the Laboratory 

(14) format 

Subgroup 2 - mass point locations {r^}. One card for each mass point 
(3E15.8) format 

Subgroup 3 - masses, Inertia matrices 

(repeat for each mass point) 

Card 1 - mass nu [(E15.8) format] 

Caxd 2 - column 1 of inertia matrix [1^] 

Card 3 - column 2 of inertia matrix 

Card 4 - column 3 of inertia matrix 

(3E15.8) format 

Group c 

(if IFLAB =0, skip to Part 5) 

Subgroup 1 - (One card) 

Card 1 - enter p (number of modes ) 

(l4) format 

Subgroup 2 

enter the p modal masses 
3 per card 
(3E15.8) format 
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Subgroup 3 

enter the p frequencies 
3 per card 
(3E15.8) format 
Subgroup 4 

enter the modal matrix - column by column 
Use 3 elements of a column per card 
(3E15.8) format 

Group d 

(if IFLAB = 0 skip to Part 5) 

Subgroup 1 

(if JUMP = 1 skip to Subgroup 2) 

enter the elastic displacements and velocities 

all {q^}, one card per mass pt [(3E15.8) format] 

all {q.}, one card per mass pt [(3E15.8) format] 

all {0^}, one card per mass pt [(3E15.8) format] 

all one card per mass pt [(3E15.8) format] 

Subgroup 2 

enter the ratio of damping to critical damping for each mode. 
3 per card 
(3E15-8) format 

PART 5 - Rigid-body displacements and velocities 
Card 1 

U] 

(3E15.8) format 

Card 2 

{h} 

(3E15.8) format 

Card 3 

M 

(3E15.8) format 
Card 4 {w X } 

(3E15.8) format 
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FART 6 - Fluids and moving mass input 
Group a - reservoir input 

(skip to Group c if no fluids axe present) 
Subgroup 1 - (one card) 

Card 1 - entry p, A 

(2E15.8) format 

Subgroup 2 - emptying reservoir 


Card 1 - 

•? 

(E15.8) format 

Cards 2-4 

[i?] 1 row per card 

(3E15.8) 

format 

Card 5 

[s?} 

(3E15.8) 

format 

Card 6 

r Ri 
{ e i3 

(3E15.8) 

format 

Card 7 


(3E15-8) 

format 

Card 8 

{u®} 

(3E15.8) 

format 

Caxd 9 


(3E15.8) 

format 

Cards 10- 

12 [<J?] 1 row per card 

(3E15.8) 

format 

Card 13 

^i* ^io* h i max’ ^i min’ 

(4E15.8) 

format 

Card l4 

; {gJ] 

(3E15.8) 

format 


Subgroup 3 - filling reservoir 
same as Subgroup 2 
Group b - input for through pipes 

(skip to Group c if no fluids or if no pipes) 


REPEAT for each pipe. 

Card 1 m.' 

i 

Card 2 {S^} 

Card 3 {Au^} 

Card 4 {A.e^} 

Card 5 } 

Card 6 { AIR } 

Cards 7-9 [4<R ] 1 row per caxd 

Cards 10-12 ■[! ] 1 row per card 


(E15.8) format 

(3E15.8) format 
(3E15.8) format 
(3E15.8) format 
(3E15.8) format 
(3E15.8) format 
(3E15.8) format 

(3E15.8) format 
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Group c - moving masses input 

(skip to part 7 if no moving masses are present) 
enter for each moving mass present 
1 per card 
(E15.8) format 


PART 7 - Counterweight input data 

same as instructions for PART 4, but for Counterweight 
Group a (one card) 


Card 1 - IPHS1, JUMP 
(214) format 

IPHS1 - jO no Phase I tape provided 
'l Phase I tape provided 
JUMP = 1 is used to set the elastic displacements 
and velocities = 0 without reading them in. 


Group b 

(If IPHS1 = 1, skip to group d) 

Subgroup 1 (one card) 

Card 1 - enter the no. mass points on the Counterweight 

(l4) format 


Subgroup 2 

mass point locations {r }. One card for each mass point 

SL 

(3E15.8) format 

Subgroup 3 - masses, inertia matrices 
(repeat for each mass point) 

Card 1 

mass m_ 
a 

(E15.8) format 

Card 2 - column 1 of inertia matrix [l_l 
a 

Card 3 - column 2 of inertia matrix 


Card 4 - column 3 of inertia matrix 
(3E15.8) format 
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Group c 
(If IFCW 
Subgroup 1 
Card 1 - 

Subgroup 2 


Subgroup 3 


Subgroup 4 


Group d 
(if IFCW 
Subgroup I 

a11 

all {qj 

all {e & } 

• 

all {9a} 
Subgroup 2 


= 0, skip to PART 8) 

- (one card) 

enter p (number of modes) 

( 14 ) format 

enter the p modal masses 
3 per card 
(3E15.8) format 

enter the p frequencies 
3 per card 
(3E15.8) format 

enter the modal matrix - column by column. 

Use 3 elements of a column per card 
(3Ei5.8) format 

= 0, skip to PART 8) 

(if JUMP = 1 skip to Subgroup 2) 

enter the elastic displacements and velocities 

, one card per mass pt (3E15.8) format 

, one ; card per mass pt (3E15.8) format 

, one card per mass pt (3E15.8) format 

, one card per mass pt (3E15.8) format 

enter the ratio of damping to critical damping for each mode. 
3 per card 
(3E15.8) format 
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PART 8 - Connecting Structure input data (skip to part 9 only if 
IFCN = 0 and KLUE = 4) 


Card 1 

i, a 



(214) format 


Card 2 

ts] 

(3E15.8) format 

Card 3 

[s] 

(3E15.8) format 

Card 4 

{6} 

(3E15.8) format 

Card 5 

{6} 

(3E15.8) format 

Card 6 

{0*} 

(3E15.8) format 

Card 7 

{**} 

(3E15.8) format 

Card 8 


(3E15.8) format 

Card 9 

} © 

(3E15.8) format 

PART 9 - 

Type 2 controls input 



Group a - (one card) 

Card 1 - enter f£ 1 
L o J 

(3E15.8) format 

Group b - pump command (PMPCOM) input (skip to Group c if there are 
no fluids present) 

Card_l - V, V[J , tfu , acc mG 

(4E15.8) format 

Group c - moving mass command (MUCOM) input for Mass 1 (j = l); 

elevator (skip to Group e if moving mass j = 1 is not 
present) 

Subgroup 1 

Card 1 - enter 3 clues, 0 or 1, denoting absence or presence of 
motion along X^, X^, and X^ axes respectively 

(312) format 

Card 2 - enter initial value of {U.} 

(3E15.8) format 

Subgroup 2 

Repeat Card 1 below for each axis present. If the mass does 
not move, skip to Group d. 

Cardi - enter acc^, vel^ 


(4E15.8) format 
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Group d - Balance mass control input (controls motion of moving 
mass j = 2) (skip to Group e if moving mass j = 2 is 
not present ) 

Card 1 - {U 2 } 

(3E15.8) format 
Card 2. - {if/ 

(3E15.8) format 

— rd 3 ~ ^1EST* K eXSE» ^DESE* T RAL» “lC BAL 
(5E15.8) format 
i 

Card 4 - {X A S } 

( 3E15 . 8 ) format 

" ^SENS/CM^ 

(3E15.8) format 

Group e - moving mass command (MIJCOM) input for masses 3, 4, ... 

(j > 2) (skip to Part 10 if there are no other moving 
masses present) 

repeat input for each moving mass as described in Group c 

PART 10 - Beam- type connecting structure, numerical integration, and 
printout data 


Group a - beam connecting structure subroutine input 


(skip to Group b if KLUE 

« 4) 

Card 1 - M c , M l , I c 

(3E15.8) format 

Card 2 - I, * J_, J T 
L’ C L 

(3E15.8) format 

Card 3 - a, b,.-AE 

( 3 EI 5 . 8 ) format 

Card 4 - El, GJ. 1 

(3E15.8) format 

°MIN 


Card 5 - V Y p , V T 



Group b - numerical integration input - see MAGIC writeup in 
Appendix B2 

Card 1 - DELTA, EPSIL, STEP 
( 3 EI 5 . 8 ) format 


* For the control systems provided with this program, the first two 
components of {Ug} must be read in as zero. 
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Card 2 - DTAU, NTIME 

(E15.8, 13) format 

DTAU is the MAGIC output time interval. In this program 
a distinction is made between the plot output interval 
and the print output interval. DTAU is used as the plot 
output interval and NTIME * DTAU is used as the print 
output interval. 

Remaining Cards - enter ERR parameters. 1 parameter for each 

integrated variable. 3 parameters per card 
(3E15.8) format 

The order of the integrated variables Y is as follows: 





l r! 
iw‘l 

iji 

?aR j 

f^ Y i 

{R\ 

w 

, t'K} 

•< 

fARj 
{**) 

K 

&Cr 

t,r «rr 


> 


rr 


U 

U 


ul 




M 


FI L T 


always present 

present if Laboratory is flexible 
present if Space Station is flexible 

present if Counterweight is flexible 
always present 

present if Laboratory is flexible 
present if Space Station is flexible 

* 

^ present if Counterweight is flexible 


(wobble damping parameters) present if 
f controls are present 


(3.5) 


l (balance mass parameters) present if 
r mass balancing is present 
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Group c - flexible problem printout data (skip if KLUE >3) 

Subgroup 1 - (one card) 

Card 1 - KLOUTL, KLOUTC 
(213) format 

KLOUTL = number of mass pts for the Laboratory at which 
the elastic displacements and velocities are 
printed each output cycle 

KLOUTC = similar to KLOUTL but for the Counterweight 
Subgroup 2 - (skip if KLOUTL = 0) 

enter the index numbers of the KLOUTL mass points 
20 per card 
(2013) format 

Subgroup 3 - (skip if KLOUTC = 0) 

enter the index numbers of the KLOUTC mass points 
20 per card 
(2013) format 

3.7 PROGRAM OUTPUT 

Examples of the program output are presented and discussed in Appendix B2. 
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APPENDIX A 

DFTATT.TO INFORMATION 
FOR THE PHASE I COMPUTER PROGRAM 


This appendix contains detailed programming information for the Phase I 
computer program. In addition, samples of input and output data are presented. 
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APPENDIX A1 
EXAMPLE OF INPUT DATA 


The example on the following pages shows input data for the 9 -m °dule 
cantilevered configuration shown in Figure Al. Module 1 is the core module 
described in Appendix F of Volume 1, and the remaining modules are identical 
and have the properties of the typical appended module described in that 
appendix. The local x axis of each appended module is oriented along the 
module, and the local z axis is oriented as indicated in the table on 
Figure Al. The last mass point of the idealized core module is fastened to 
the ground and therefore does not enter the problem. Mass points 2 and 5 
have zero mass and were artifically introduced to serve as connection points. 
20 modes of Module 1, and 10 modes of all other modules are used in the syn- 
thesis. Mass-loaded modes are used for Module 1 (see Section 2.2). As the 
T-matix data indicates, the x-direction deflections in the input modes of 
masses 7 and 8 are not present. These deflections were constrained to be 
zero in the core-module idealization. All other mass points have six degrees 
of freedom. 

The results of this problem include printout of all input substructure 
modes (see Appendix A2). 
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Fig. A1 Configuration Used to Illustrate Input Data for Phase I Program 
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MASS POINT LOCATIONS 
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APPENDIX 2 

EXAMPLE OF PROGRAM OUTPUT 

The computer printout of the problem described in Appendix A1 appears 
on the following pages. The illustration of Appendix A1 indicates that 100 
frequencies are to be computed; however, the input was revised so that only 
20 frequencies were computed. 
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APPENDIX A3 

BRIEF DESCRIPTION OF EACH SUBROUTINE 


CONFIG - 
LOGLO - 

AIJMOD - 
MASS 

MODE 

TDEF 

FILE 

JOINT - 

FORMB - 

BTMB 

RTDS 

MXI 


INCH 

DIAG 

OFF 

LATE 

FREE 

El GOUT - 
MPRD 
TPRD 
MCFY 
L0C810 - 
L0C811 - 
ZERO 
ADD 


Checks to see that the program limitations are met 

(For each module) Reads the mass point locations in local coordinates 

and transforms them to global coordinates 

(For each module) Forms the [A. .] matrix 

ij 

(For each module) Reads the physical mass matrix in local coordinates 
and transforms to global coordinates 

(For each module) Sets up T-matrix constraints, reads mode shapes in 
local coordinates and transforms to global coordinates 
Sets up each T-matrix 

Rewrites the file containing the physical mass matrix in global 
coordinates for use in later processing 

Creats a file of mode shape partitions and one of [A^j ] matrix 
partitions at junction points for use in forming the B-matrix 
Forms the B-matrix 
Forms B T MB 

Used by BTMB to form R T DS where D is a partition of the block diagonal 

mass matrix and R and S are partitions of the B-matrix 

Reads in modal mass and eigenvalue matrix, stores B^MB partitions into 

M^. matrix and writes it out. MXI then stores the stiffness matrix in 

the proper form (compressed) for FUTILE. MXI will call CHEAP to perform 

free-free conversion of problem if needed 

Used by CHEAP to invert a symmetric matrix 

Used by MXI to set up partitions of M^ matrix lying along diagonal 

Used by MXI to set up partitions of M^ matrix lying below the diagonal 

Solves the eigenvalue problem and orthonormali ze s the vectors 
Back transformation of vectors for free-free systems 

Transform the eigenvectors back to the physical global coordinate system 
Multiply two matrices 

Multiply the transpose of a matrix by another matrix > 

Copy one matrix into another 
Locate member of file 10 - B partition 
Locate member of file 11 - global mass partition 

Zeroes out a matrix 
Obtains the sum of two matrices 
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SEQUENCE OF OPERATIONS IN SUBROUTINE LATE 

Subroutine LATE requires additional discussion. This subroutine solves 
the coupled-structure eigenvalue problem which is in the following form: 

AX = \BX where A is symmetric 

and B is symmetric, positive definite 


CALL RJTILE to find the Cholesky Decomposition of B 

T 

B = L L 


-1 -T 

CALL DAGGER repeatedly to form L A L 

the eigenvalue problem to be solved is 

( L -1 A L" T )(L T X) = \ (L T X) 

-1 -T 

CALL SWITCH to store the matrix L A L in a form acceptable to SYMEIG 

CALL SYMEIG to solve the eigenvalue problem giving the X’s 

CALL SYMVEC (An entry point of SYMEIG) to give the eigenvectors Y = L T X 

CALL ORTHON to orthonormalize the set of eigenvectors Y 
CALL TRIEQ to find X given Y = L T X 


For this program 

A is the matrix, and 

B is the k^ matrix except for nonsingular systems (those with rigid-body 
modes) where 


A = ([** 22 -^; " ^12-^ ^ M ll^ and 
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APPENDIX A5 

INPUT, OUTPUT, AND I NTE MEDIATE FILES 

FILES 1, 2, 3, 4 , 7 , 8, 9, 10 - intermediate files 

These files are used to temporarily hold matrix partitions to be used 
later in the program. All files are unformatted. Details are indicated in 
the program flow charts. 

FILE 11 - output file 

File 11 is the output file to the Phase II program corresponding to 
run number 1 (usually the Laboratory run). Write statements are in CONFIG, 
LOGLO, FILE, and EIGOUT. 

FILE 12 - output file 

File 12 is the output file to the Phase II program corresponding to 
run number 2 (usually the Counterweight run). The write statements are the 
same write statements as those used in FILE 11. 
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PROGRAM FLOW CHARTS 


MODID 

NMO 

NMPO 

IXYZ 

G 

NIX 

NUMEIG 

NUMVAL 

NZERO 

XYZ 

XYZBAR 

MPN 

NMP 


Symbols Used in the Flow Charts 

The array used to denote absence or presence of each module. 

The array used to keep the number of modes for each module 
The array used to keep the number of mass points for each module 
The array used to store mass point numbers of the attachment points 
for each attached module 
Constant G 

Error return indicator from INCH and FUTILE 
Number of eigenvectors to be found by SYMVEC 
Number of eigenvalues to be found by SYMEIG 
Number of rigid body modes 

Array containing X^, X^, Xg, and mass point locations when read in, 
then later holding mass point locations in absolute global coordinates 
Array containing mass point locations in relative global coordinates 
when computed 

Current mass point number for a module 
Number of mass points for a module 



MAIN PROGRAM 
























































SUBROUTINE TDEF 















SUBROUTINE FORM B 























Rewind 8 
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t/5 






















REWIND 8 



(6did( if) 
ffiSo(i] 





















STORE THE CO 
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SUBROUTINE CHEAP 
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SUBROUT IKE FREE 
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SUBROUTINE EIGOUT 



ACCUMULATE MULTIPLICA- NORMALIZE TRANSFORMED VECTOR TO A 

TION BY ROW IJ OF B UNIT MODAL MASS 
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APPENDIX A7 

FORTRAN SOURCE DECK LISTING 


The fortran source deck listing appears on the following pages, 
listing shows the LRC CDC version of the program. 


This 



PROGRAM SP ACE( INPU T= 20 1 ,QU TP'JT = 201 , T A PE5= l NPUT , TAPE6=3JTPJT , 100000 

TAPE 1=401, TAP E2= 401 , TAPE 3=401 , TAPE4 = 40l , 200000 

TAPE8,T APE9 = 401, TAPE 10, TAPE 11 = 401 , TAPE 12 = 201 1 300000 
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180 KEY = 0 MXI07740 167800000 

GO TO 100 MXI07750 167900000 

190 RETURN MXI07760 168000000 

END MXI07770 168100000 
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GO TO 145 18X1 08720 175900000 

170 RETURN MXI08730 176000000 

END MX I 08 740 176100000 



SUBROUTINE PUT l L E ( A , N, Nl X » MXI08750 176200000 
DIMENSION AIL) MXI08760 176300000 
DOUBLE PRECISION SUM MXI08770 176400000 
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APPENDIX B 

DETAILED PROGRAMMING INFORMATION FOR 
THE PHASE II COMPUTER PROGRAM 


This appendix contains detailed programming information for the Phase II 
computer program. In addition, samples of input and output data are presented 
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APPENDIX KL 

EXAMPLES OF INPUT DATA 


The following two examples illustrate the method of setting up the input 
data deck. The configuration for both examples is that discussed in Section 6 
of Vol ume I. Computer output for these examples is presented in Appendix B2. 

Example 1 

This example shows input data for a 40 second rigid body run with plotting. 
This run has fluids present, no deployment, no controls, no supplementary 
forces, no moving masses and no structural loads computed. Although flexible- 
body data is present, the constraint option has been used to rigidize the 
entire Space Station. The JUMP clue to set all {q.^} , {q^}>{0^} and {9^} to zero 
is used. Fixed- step numerical integration with a step size of .05 is indicated 
together with a plotting output interval of .09997436 so that plotting will 
occur the first time a multiple of .05 exceeds each multiple of .0997436. 

The time-history printer output interval is 2* .0997436. 



final time, pilot scale 
plot clue, fluid clue, no. pipes 
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DELTA, EPSIL, STEP 
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Example 2 

This example shows input data for a 10 second run with with KLUE = 3 
(flexible Laboratory and rigid Counterweight). The quiescent-state initial 
conditions are used for the deformations (see Section 3.4.10). The complete 
quiescent-state initial conditions are presented in the output shown in 
Appendix B2. Plotting will be done, and moving masses 1 and 2 are present 
indicating elevator motion and mass balancing. Controls are present, and 
loads will be printed. The integration step-size is fized. 



final time, plot scale 
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APPENDIX B2 


EXAMPLES OF PROGRAM OUTPUT 

This appendix shows the output for the two problems described in 
Appendix Bl. The key to reading the titles of the time-history print out 
is given at time = 0 in these examples . 

Example 1 

The first portion of the output consists of the input data. The input 
modes were truncated; however, they are shown in their entirety in Appendix 
I of Volume I. 

The second portion of the print out is the time history. Auxiliary 
results presented in this portion include the location of the center of mass 
and the total-system angular momentum. 
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Example 2 

The time-history printout for this problem, where the vehicle is flex- 
ible, includes the connecting-structure deformations {6} and {0*} as well 
as the main variables. Elastic deformations are also printed at each mass 
point requested. In addition, internal structural loads are printed at the 
requested points. As discussed in Section 4.5 of Volume I, the internal loads 
are inaccurate at t =0. These loads have been deleted from the. printout at 
t = 0; however, this modification was done after making the run on the 
example presented. After completing the time history, a tabulation is 
printed of the maximum and minimum (or largest negative) elastic motions 
that occurred during the run at every mass point. Also, the maximum and 
minimum internal structural loads are printed at only those points where 
loads were requested in the input data. The maxima and minima for both 
deformations and loads are obtained by examining the data at every print- 
out time point and not at every numerical-integration time point. 
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APPENDIX B3 

Description of Numerical- Integration Subroutine 
(MAGIC) Parameters 

The numerical-integration subroutine provides the user with both a 
fixed-step and variable-step numerical-integration option. Both options 
use a Runge-Kutta method. For details of this subroutine, the user is 
referred to the MAGIC documentation manual.* 

The parameters needed to control this integration process are DELTA, 

EPSIL, STEP, DTAU, and ERR. 

B3.1 Fixed- Step Runs 

For fixed-step runs set DELTA = 0 and set STEP to the time step size of 
each numerical-integration cycle. A guideline for -choosing STEP is to set 
it to l/20th of the highest frequency present among the integrated variables. 
DTAU is the time interval between both the print and plot output cycles. 

EPSIL and ERR are not used in fixed-step runs. 

B3.2 Variable -Step Runs 

In variable-step runs, the subroutine tests the accuracy of the integra- 
tion for each variable and uses this information to set the numerical-integra- 
tion time step size. As explained in detail below, DELTA and EPSIL are small 
positive numbers used to test the integration accuracy by comparing a precision 
parameter with a percentage of each variable. When the absolute value of a 
variable is very small an ERR parameter (one is supplied for each integrated 
variable) is used for testing accuracy. One guideline for selecting DELTA and 
EPSIL is to choose DELTA = l/3 EPSIL. The approximate tests started below 
describe the usage of EPSIL, DELTA, and ERR and should be used to determine 
their values^. MAGIC initially selects a step size of . 0078125 . 


* Ball, D. , Stoodley, G. , and Angelillo, R. , "A Fortran Program for the 
Numerical-Integration of Ordinary Differential Equations," Grumman Aerospace 
Corporation, Data Systems Department Report, CR-61-2, October 1969 . 
t The precise test is more complicated and is stated in the previously 
referenced MAGIC documentation manual. 
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Let y be a typical element of the vector {y} to be integrated. Then 
y dev is a precision parameter determined by subtracting the values of y 
obtained by using two somewhat different numerical-integration processes. 

The tests are as follows: 

Dominant test when y is not small: 

If|y dev | > EPSIL*|y| for any integrated variable, the step 

size is halved and the integration is 
repeated 

If|y dev| < DELTA*! y | for each integrated variable, the 

integration is not repeated and the step 
size doubles 

If | DELTA*! y | <|y dev|< EPSIL*]y| integration continues with the current 

step size 

Dominant test when y is small (of order ERR): 

If|y dev|> ERR for any integrated variable , the step 

size is halved and the integration is 
repeated 

If|y dev | <.01 ERR for each integrated variable, the integration 

is not repeated and the step size doubles 
If .01 ERR ^|y dev|^ ERR integration continues with the current 

step size 

Discontinuities in the analysis (e.g., pulse-type forces or an elevator 
moving with constant acceleration) may cause the step size to halve indefinitely. 
To avoid this difficulty, the user must input a value of STEP which is used as 
the minimum time interval. STEP has no other use in variable step runs. 

DTAU is the time interval between print output cycles. 

B3.3 Limitations of variable-step Option 

While the variable time-step option often reduces the computer running 
time considerably, the following problem limits the use of this option. 
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Whenever the step size is halved all computations needed to obtain the 
derivatives are repeated at a value of time which is less than the previous 
valued used for these computations. Several control-system subroutines 
use the value of the variables at the previous time step; however the stored 
value now corresponds to a future time step. This discrepancy leads to 
errors when these subroutines are used. For this reason the following prob- 
lems should not be run using the variable-step option: 

Moving masses and mass balancing 
Fluid motion 

Deployment or retraction maneuvers 
Spin up or spin down 
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APPENDIX B4 

BRIEF DESCRIPTION OF EACH SUBROUTINE* 


ADD 

ATTCTR 

BMCON 

BMCTR 

CHASE 

COMMIN 

CONSTR 

CONTRL 

CW 

CWCOM 

DOTP 

FIRST 

FUTILE 

GAMMA 

INCH 

LAB 


LAST 

LOADIN 

LOOK 

LUMP 

MCPY 

MINCY 


MINMAX 


addition of 3 X N matrices 
attitude control subroutine 
balance mass control interface routine 
balance mass control subroutine 

solution of equations for symmetric positive definite 
system 

major subroutine for data input and initialization of 
common areas 

connecting structure subroutine for circular beam 
interface program for controls subroutines 
calculates the contributions to [ A ] and { V } 
from the Counterweight 
Counterweight command 
dot product of vectors 

computes the realignment of the initial conditions 
Cholesky decomposition of symmetric positive definite 
matrix 

compute gamma functions [see Appendix B, Volume I] 
inversion of a real symmetric positive definite matrix 
calculates the portions of [ A ] and { V } corres- 
ponding to the Laboratory, not including fluids and 
moving masses 

compute connecting structure deformations 
computations of structural loads 
look-up subroutine to give the values of k denoting 
the masses m^. which influence moving mass 

calculates the contributions to [ A ] and { V } 
from the moving masses 
matrix copy subroutine 

sets up and performs integration step, prints time - 
history output and plotting output and computes the 
transformations to physical coordinates 
computation of minimum and maximum displacements, 
velocities and accelerations 


* For the definition of the analytical symbols, see Section 4 
of Volume I. 
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MPKD 

MTRAN 

MUCX)M 

MULT 

PIF 

PIPE ' 

PL03AX 

PMPCOM 

POSCOM 

POSCTR 

FTLP 

RESVR 

SCALR 

SETUP (MAGIC) 
SPNCOM 
SPNCTR 
SSIN 

SUBTR 

SUPPLM 

TPRD 

VELCOM 

VELCTR 

WBLCTR 

XADD 

XMCPY 

XMPRD 

XTPRD 

ZCW 

ZLAB 


general matrix molt ipli cation 
multiplication by transpose of a 3 x 3 matrix 
moving mass command, subroutine 
multiplication by a 3 x 3 matrix 
computes n_ function [see Appendix B, Volume I] 
calculates the contributions to [ A ] and { V } 
from the pipes 
Calcomp plotting program 
pump command subroutine 
position command subroutine 
position control subroutine 
T 

compute [P] [L] [P] where [P] is a 3 x 3 matrix and 

L is the 3x1 diagonal of a 3 x 3 diagonal matrix 
calculates the contribution to [ A ] and { y } from 
the reservoir 

multiplication of a matrix by a scalar 
Runge-Kutta integration program 
spin command subroutine 
spin control subroutine 

computes S-function or S" 1 function [see Appendix B, 
Volume 1] 

matrix subtraction subroutine 

supplementary forces subroutine 

general matrix-transpose multiplication 

velocity command subroutine 

velocity control subroutine 

wobble control subroutine 

general matrix addition 

matrix stacking subroutine 

general matrix multiplication 

general multiplication by transpose of a matrix 

calculates portions of [A r ] and { y R ] 

calculates portions of [A r ] and { V R 3 
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APPENDIX B5 

REPLACEABLE SUBROUTINES 

The following subroutines may be replaced by the user in order to run 
different physical problems: 

CONSTR - computes structural loads for beam-type connecting structure 

CONTRL - main controls interface routine; calls and computes data for CWCOM, 

SFNCOM, ATTCTR, SPNCTR, and WBLCTR 

CWCOM - commands position of Counterweight 

SPNCOM - commands Space -Station spin speed 

MUCOM - commands motion of moving masses (j^.'s) 

0 

PMPCOM - commands velocity of fluid being pumped 
POSCOM - position command policy 
VELCOM - velocity command policy 
ATTCTR - controls Space -St at ion attitude 
SPNCTR - controls Space-Station spin speed 
WBLCTR - controls Space-Station wobble 

BMCON - interface for BMCTR; calls and computes data for BMCTR 
EMCTR - controls motion of balance mass to balance mass shifts along long 
axis of Space-Station 
POSCTR - position control policy 
VELCTR - velocity control policy 

SUPPLM - contains supplementary external loads as a function of time (these 

are set to zero in the subroutine supplied) 

LOOK - selects structural masses (m^'s) that influence motion of moving 

masses ( p. . ’s). In the version of LOOK supplied with the program 

k = 2, 8 for all (i . ' s and the k's remain constant with time 

J 

CONTRL and BMCON are interface subroutines for the major control systems 
(i.e., the control systems required to control the motion of the Space Station). 
Figure B1 illustrates the hierarchy of the control- system subroutines. 

The read and print statements for the replaceable subroutines are 
located within these subroutines. 

The input and output symbols for all replaceable subroutines are 
tabulated on the following pages. The analytical symbols used in this 
chart are defined in Section 3.2. 




Fig. B1 Hierarchy of Control Subroutines 
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CONSTR - Called from ZLAB and ZCW 

INPUT 


Main Program Analytical Symbols 

CONSTR Symbols 

Comments 


EL 

( average axial load compo- 

{P ^ } 3 

P 

\ nent in Connecting St rue- 


[ ture 

{$*} 

PHI ST 


{$*} 

PHI STD 


{6} 

DEL 


{6} 

DELD 


t 

T 



IN 

Card Input Unit 


IOUT 

Printer Output Unit 

OUTPUT 



Main Program Analytical Symbols 

CONSTR Symbols 

Comment s 
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CONTRL - Called from ZIAB and ZCW 

INPUT 


Main Program Analytical Symbols 

CONTRL Symbols 

Comments 


N 

Number of mass points on the vehicle 


IFCON 

(Set to 1 if controls are present 
lue | Set to 0 if controls are not present 


LABCW 

( Set to 1 if called frcm ZLAB 
Llue 1 Set to -1 if called from ZCW 

t 

T 


{u) X } 

W X 


all fr. } 

RI 

{r_^} begins at RI (3 • i - 2) 

all {C^} ‘ 

TH 


all {0.} 

THD 


(R] 

RDDOT 


,*X, 

(Cl) } 

WXDOT 


all {c*.} 

XDD 

begins at XDD (1, i) 

ail fe i } 

XDD 

ffi } begins at XDD (4, i) 

M c 

EMC 

mass of Counterweight (see Note 6) 

M 

TOT 

mass of Space Station 

“l 

EML 

mass of Laboratory (see Note 6) 


GAM 


fy] 

GAMD 

— 


EL 


{i 0 } 

ELD 


c'V 

ELDD 

' 


IN 

Card Input Unit 

r 

IOUT 

Printer Output Unit 

0 


1 1 

t* 
t r» 

i i 

DL 

> See Note 6 

i — i 

w 

o 

1—1 

DC 

f 

Pr.il 

IL 

Inertia matrix of Laboratory 

1 — 1 
H 
Q 
i i 

IC 

Inertia matrix of Counterweight 


YPRED 

Integrated controls variable (after t = 0) 


Continued 
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CONTRL (Continued) 


OUTPUT 


1 Main Program Analytical Symbols 

CONTRL Symbols 

Comments 1 

all Cf 1 } (or (f^l ) 
c c 

A 

(f^j (orj{f^} )begins at A(3i -2) (see Note 1) 

all CtJ] (or Ct*}) 

A 

(t^} (or(f^} )begins at A(3i +l)(see Note 1) 

r'X^S, 

} 

DO 

Output when entered via EARLY 

i 

s 

IS 



NEQS 

The number of integrated controls variables 



(at T = 0) 


YICS 

Initial conditions of integrated controls 



variables (at T = 0) 


YDOT 

Time -derivative of integrated controls 



variables 


Special Considerations 


(1) The forces and torques exerted by the control actuators are computed for the: 

Laboratory when LABCW = 1 
Counterweight when LABCW = -1 

(2) Entry point CONTRD is used to input the controls input data 

(3) Entry point EARLY is used to compute the value of 



(4) Entry point LLDLDD is used to call CWCOM and then return 

(5) The variables currently being integrated are 4 variables from WBLCTR 
(if controls are present) and 3 variables from BMCTR (if mass balancing 
is present) 

( 6 ) See Equations (A29) and (A 36 ) of Appendix C, Volume I 
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CWCOM - Called from CONTRL 


INPUT 


Control - System 
Analytical Symbols 

CWCOM Symbols 

Comment s 

t 

TIME 


^ 03 u 

XB3CU 

J 

Connecting Structure length 
update 

t I03n 

TXB3CU 

Update time 

9/C C .. . 

103 max 

ACCB3C 

Magnitude of acceleration of 
t 0 3 

vel 

i !03 max 

VEIMBC 

Maximum velocity of 

^3 

XB3C 

Input at t = 0 

OUTPUT 

Control - System 
Analytical Symbols 

CWCOM Symbols 

Comments 


DDXB3C 


' £ 03 

DXB3C 


£ 03 

XB3C 

Jndeformed length of connect- 
ing Structure 
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SPNCOM - Called from CONTKL 


INPUT 


Control - System 
Analytical Symbols 

SPNCOM Symbols 

Comments 

time 

TIME 


CD 

1CU 

01CU 



T01CU 


0)1 CU 



acC 0)l C 

ACC01C 

At time = 0 

f 

FJET 

| These variables are 

JET 

A. 

(X 

I in the argument list 


XAJ 

l but are not currently 

[Mu] 

WM11 

' used. 

“ic 

01C 

At time = 0 


OUTPUT 


Control - System 
Analytical Symbols 


SPNCOM Symbols 


Comments 
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MUCOM - Called from LUMP 


INPUT 


Main Program Analytical Symbols 

CONSTR Symbols 

Comments 

t 

T 


j 

L 

Moving mass number 


IN 

Card input unit 


IOUT 

Printer output unit 


OUTPUT 


Main Program Analytical Symbols 


CONSTR Symbols 


Comments 
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FMPCOM - Called from MAIN 


INPUT 



OUTPUT 
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POSCOM - Called frcm MUCQM, CWCOM 


INPUT 


Control - system 
Analytical Symbols 

POSCOM Symbols 

Comments 

time 

TIME 


poscu 

POSCU 

Position command update 

timcu 

TIMCU 

Time command update 

acc 

mag 

ACCMAG 

Magnitude of acceleration command 

vel 

max 

VELMAX 

Maximum velicity command 

OUTPUT 


Control - system 
Analytical Symbols 

POSCOM Symbols 

Canments 

out 

OUT 


timeO 

TIMEO 

\ 

timel 

TIME1 


time2 

T3ME2 


time3 

TIME3 


accOl 

ACC01 

Intermediate output to each subroutine 

accl2 

ACC12 

'that calls POSCOM 

acc23 

ACC23 

[and returned each time cycle 

poscO 

POSCO 


poscl 

P0SC1 


posc2 

P0SC2 


veld 

VELC1 


velc2. 

VELC2 


acc 

cmn 

ACC CMN 

Acceleration command 

vel 

cmn 

VELCMN 

Velocity command 


POSCMN 


Position command 










VELCOM - Called from SFNCOM, EMPCOM 


INPUT 


Control-system 
Analytical Symbols 

VELCOM Symbols 

Comments 

time 

TIME 


velcu 

VELCU 

Velocity command update 

timcu 

TIMCU 

Time command update 

acc 

ACCMAG 

Magnitude of acceleration command 

mag 
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ATTCTR - Called from CONTRL 


INPUT 


Control- system 
Analytical Symbols 

ATTCTR Symbols 

Comments 

(0} 

T 


C©3 

DT 


e e 

DB 

ETDB 


• 

MAX 

DETMAX 


A 

t-x 

XAJ 


1 1 

H 

i i 

WM11 

Inertia matrix of Space Station 

f jet 

FJET 


'"iCMAG 

OICMAG 



OUTPUT 


Control- system 
Analytical Symbols 

ATTCTR Symbols 

Comments 

A. 
f J 






FAJ 

CONTRL currently 
expects values for: 

FAJ 

(2, 

1)* 





FAJ 

(2, 

2) 




FAJ 

(2, 

3) 




FAJ 

(1, 

b) 




FAJ 

(2, 

5) 




FAJ 

(2, 

6) 




FAJ 

(2, 

7) 




FAJ 

(1, 

8) 




FAJ 

(2, 

9) 




FAJ 

(2, 

10) 




FAJ 

(2, 

11) 




FAJ 

(1, 

12) 




FAJ 

(2, 

13) 




FAJ 

(2, 

lb) 




FAJ 

(2, 

15) 




FAJ 

(1, 

16 ) 


* The first index designates the axis and the second designates the jet number 
(1 through 16) indicated in Figure 3.8. FAJ(I,J)=0 when jet is off. Otherwise 
FAJ(l,J)=+ Fj e £ depending whether the jet is directed along the positive or 

negative jth axis. 











SPNCTR - 


INPUT 

Control - System 
Analytical Symbols 

SPNCTR Symbols 

°lc 

OIC 

(»} 

0 

e .. 

E01DB 

uxLDB 

f 

FJET 

JET 


“icMAG 

01CMAG 

OUTPUT 


Control-System 

Analytical 

SPNCTR Symbols 

fAj 

FAJ 
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from CONTRL 


Comments 


Comments 


See Table for ATTCTR on previous 
page 
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WBLCTR - Called frcm CONTRL 


INPUT 


Control - system 
Analytical Symbols 

WBLCTR Symbols 

Comments 

C to} 

0 


M 

WM11 


P 

BET 

equivalenced to YPRED (2) 

P 

DBET 

equivalenced to YPRED (l) 

*t 

KT 


t 

TQLIM 


qLIM 


■S 

WKDBET 


Co)} 

DO 


time 

TIME 


mom CMG 

WMOMCM 


T 

WIBGIM 


PGIM 


err KT 

ERRKT 


err KTI 

ERRKTI 

equivalenced to YPRED ( 3 ) 

• • 
P 

DDBET 


“ICMAG 

OICMAG 


err LL 

ERRLL 


err LLI 

ERRLLI 

equivalenced to YPRED (4) 

i 

IS 


s 




YICS 

initial conditions of integrated 
controls variables 


YPRED 

integrated controls variables (after T=0) 


NEQS 

the number of integrated controls 
variables 


Continued 





B-64 


WBLCTR (Continued) 


OUTPUT 


Control - system 
Analytical Symbols 


WBLCTR Symbols 


YDOT 


Comments 


time-derivative of integrated 
controls variables (f3, *0, err 


err u ) 


KT 


9 


tAo 

n CMG 


TAJ 

NCMG 
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34C0N - Called frcm LUMP 


INPUT 


Main program 
Analytical Symbols 

BMC ON Symbols 

— > 

Ccmments 


IN 

card input unit J 


IOUT 

printer output unit 

all CR ± } 

RI 


all C^} 

RID 

1 

yi 

all {uj } 

WXI 

1 

{/} 

W X 

1 

1 

all {r^} 

RID 


all C0 ± J 

TH 


i s 

t/ s ) 

IS 


0 


£R} 

R 



YICS 

initial values of integrated 
controls variables 


YPRED 

gives values [Ug}^ = (5); 

£U 2 } 3 = YPRED (6) for I>0 

t 

T 


CV 

UJ 


cy 

UJD 


t^j 

UJDD 


i 



r *x, s 

DO 


U) 

OICMAG 


1CMAG ■ 



00 

RDDOT 

value computed at previous program cycle 

{^ X } 

WXDOT 

value computed at previous program cycle 


Continued 


Q 
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EMCON (Continued) 


- 

Main - program 
Analytical Symbols 

BMCON Symbols 

Comments 


ODD 

value computed at previous program cycle 

{zcm} 

ZCM 

(after T=0) value computed at previous 
program cycle (Eq. (A27b), App. C, VI) 

M 

c 

1MC 

Mass of Counterweight 

M 

TOT 

Mass of Space Station 

all 

U 


Tt (T1 ) 

PIETA 


"( Y ) 

PIGAM 


fv 

EL 


tv 

ELD 



OUTPUT 


Main - program 
Analytical Symbols 

BMCON Symbols 

Comments 

{u 2 } 

UJ 


cy 

UJD 


cy 

UJDD 



YDOT 

YDOT (5) set to [U g } ; YDOT (6) 


NEQS 

Set to {U.y set to 2 at T=0 only 


REMINDER: 

UJ, UJD, UJDD contain the values of (U^}, {u^}, £ U } upon entrance but are 

expected to contain the values of (U 0 } 5 fu_], fu 0 3 

2 ’ 2 ’ w 2 upon exit. 
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B4CTR - Called from EMCON 

INPUT 


Control- system 
Analytical Symbols * 

BMCTR Sym 

Comments 


DXB 


li} 

DO 


M 

WM y 


(\> “ 2 ) 

WMI 


i 



y S 

[r(i x )] 

GDO 


i 



[r(uo x s )] 

GO 


{x (1) 3 

DXI 

first column 

(x (2) ) 

DXI 

second column 

“lEST 

WMLEST 


K eXSE 

WKEXSE 


^ESE 

WKDESE 


T 

TAUBAL 


BAL 



(1) 

1CMAG 

OICMAG 


j- CMND i 
1 SENS/CM J 

XSCM 


[Y sens ) 3 

DDXSNS 


tx (1) } 

DDXI 

first column 

“ICBAL 

01CBAL 


t 

T 

first element 


OUTPUT 


Control- system 
Analytical Symbols 

BMCTR Sym 

Comments 

fx (2) } 

DDXI 

second column 


* See Appendix H13 of Volime I. 
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POSCTR - Called fran ATTCTR for each axis 


INPUT 


Control- system 
Analytical Symbols 

POSCTR Symbols 

Comments 

e 

E 

Position error 

e^T-. 

EDB 

Position error dead band when 

DB 



e 

EDOT 

e = 0 

®MAX 

EDO IMX 


acc 

ACCS 

Nonzero acceleration magnitude used 

s 

to establish switching curves 


OUTPUT 


Control- system 
Analytical Symbols 

POSCTR Symbols 

Comments 

acc NOW 

ACCNOW 

Multiplier of desired acceleration 
( +1. , 0. , or -1. ) 















B-69 


VELCTR - Called from EMCTR, SPNCTR 

INPUT 


Control System 
Analytical Symbols 

VELCTR Symbols 

Comments 

• 

e 

EDOT 

Velocity error 

®DB 

EDOTDB 

Velocity error dead band 

OUTPUT 

Control System 
Analytical Symbols 

VELCTR Symbols 

Comments 

acc NOW 

ACCNOW 

Multiplier of desired 
acceleration (+1., 0., or -1.) 












SUPFLM - Called from ZLAB and ZCW 


INPUT 


Main Program 
Analytical Symbols 

SUPFLM Symbols 

Comments 


N 

number of mass points on the vehicle 


K 

, , . r 1 if called from ZLAB 

-1 if called from ZCW 


IF 

11 if supplementary forces 

clue set to ) are to be computed 

jO if supplementary forces 
( are not to be computed 


OUTPUT 


Main Program 
Analytical Symbols 

SUPPIM Symbols 

Comments 

all {fj} (or {?*}) 
all (Tg] (or {t*}) 

A 

A 

{f*} (or begins at A (3i-2) 

( see Note) 

(or t T ^}) begins at A (3i+l) 
( see Note) 


NOTE: 

Supplementary forces and torques are computed for the 

Laboratory when K = 1 
Counterweight when K = -1 








LOOK - 


INPUT 


Main program 
Analytical Symbols 


LOOK Symbols 


j 


L 


OUTPUT 


Main program 


Analytical Symbols 

LOOK Symbols 


M 

k 

K 
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.ed frcm LUMP 


Comments 


moving mass no. 
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APPENDIX B6 

INPUT. OUTPUT. AND INTERMEDIATE FILES 
FILE 1 - intermediate file 

File number 1 is used for intermediate data to be used by the plotting 
program. Write statements are in MAIN and MINCY. Read statements are in 
FL03AX. 

FILE 11 - input file 

Contains Phase I input for the Laboratory. READ statements are in 
COMMIN. 

FILE 12 - input file 

Contains Phase I input for the Counterweight. READ statements are in 


COMMIN. 



Fortran 


IFLAG 

SUM 

DELT 
I PLOT 
IFPLT 
T 

IFLUID 

NPIPE 

NMLM 

NCYCLE 

FT 

KLUE 

LABCW 

YICS 

YDOT 

NTIME 

DTAU 

MAGOUT 

IFDEPL 

IFCON 
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APPENDIX B7 
FLOW CHARTS 

Symbols Used In Flow Charts* 

set to 1 when emptying reservoir becomes empty or 
when filling reservoir becomes fall 

dt where V is fluid velocity in pipe 

t 

o 

time increment 

set to 1 if plotting is present for any run on the job 
set to 1 if plotting is present for the current run 
Time 

fluid clue: set to 1 if fluids are present, set to 0 if not 

number of pipes present 

number of moving lumped masses 

number of output cycles 

final time 

- problem clue (see Section 3*3> Part 1, Group c) 

argument of CONTROL subroutine indicating whether subroutine 
was called from ZLAB (LABCW=l) or from ZCW(LABCW— l) 

- vector of initial conditions used in MAGIC integration 
vector of derivatives to be integrated by MAGIC 

| see Section 3.3j Part 10, Group 6 

Magic output clue. Magic sets MAG0UT=2 to indicate an output 
cycle 

deployment clue: set to 1 if deployment is present, set to 0 if not 

controls clue: set to 1 if controls are present, set to 0 if not 


*For analytical symbols see Volume 1, Section 4. 



MAIN PROGRAM 
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MAIN PROGRAM 























m. ii}> uiAi} 
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LAB 




INITIALIZE THE VALUES OF f 

K? *> 

J P 

i 













|l| by 
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START 


COMPUTE f X* ] , [ 

Cs 1 ! 


t 

add intom 

AND ADD j INTO \v\ 

1 



t 

COMPUTE fxH , f 
INCREMENT „ 

by 

AND K BJf 

3 

/ 

™ D r^i , [Xj] 
ADD { J INTO 

I 

, [if] INTO [ A 1 

l v J 


RETURN 








START 



C(«8) J J , $(f»2) J f 





START 


COMPUTE { ZCM J , jP-} 



CALL CONSTR TO COMPUTE 
CONNECTING STRUCTURE 
FORCES, TORQUES 


CALL CONTRL (WITH LABCW = l) TO ADD CONTROLS FORCES, 
TORQUES TO £ V } . CALL SUPPIM TO ADD SUPPLEMENTARY 
FORCES, TORQUES TO f V }. IF LABORATORY LOADS ARE TO BE 
COMPUTED: CALL LABLD 


COMPUTE PORTIONS OF 

Pj3- c/n T t<n Cat tti 

AND PORTIONS OF 

o?j T Cal { 7} 


RETURN 
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LAST 












MQCOM 





SET CONTROIS FORCES AND 
TORQUES VECTOR (A) to 0 


















bmcon 
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MINMAX 
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L0ADI1J 



sawn 
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APPENDIX B8 

FORTRAN SOURCE DECK LISTING 


The CDC LRC version of the program appears on the following pages. 




**** YOUR sequenced program follows **** 
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